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ABSTRACT 


This report contains six papers presented by the Lincoln Laboratory Air Traffic Control 
Systems Group at the American Institute of Aeronautics & Astronautics (AIAA) Guidance, 
Navigation and Control (GNQ conference on 6-9 August 2001 in Montreal, Canada. The work 
reported was sponsored by the NASA Advanced Air Transportation Technologies (AATT) 
program and the FAA Free Flight Phase 1 (FFP1) program. The papers are based on studies 
completed at Lincoln Laboratory in collaboration with staff at NASA Ames Research Center. 

These papers were presented in die Air Traffic Automation Session of the conference and 
fall into three major areas: Traffic Analysis & Benefits Studies, Weather/Automation Integration 
and Surface Surveillance. In the first area, a paper by Andrews & Robinson presents an analysis 
of the efficiency of runway operations at Dallas/Ft. Worth using a tool called PARO, and a paper 
by Welch, Andrews & Robinson presents delay benefit results for the Final Approach Spacing 
Tool (FAST). In the second area, a paper by Campbell, et al describes a new weather 
distribution system for the Center/TRACON Automation System (CTAS) that allows ingestion 
of multiple weather sources, and a paper by Vandevenne, Lloyd & Hogaboom describes the use 
of the NOAA Eta model as a backup wind data source for CTAS. Also in this area, a paper by 
Murphy & Campbell presents initial steps towards integrating weather impacted routes into 
FAST. In the third area, a paper by Welch, Bussolari and Atkins presents an initial operational 
concept for using surface surveillance to reduce taxi delays. 
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RADAR-BASED ANALYSIS OF THE EFFICIENCY OF RUNWAY USE* 


John W. Andrews 
MIT Lincoln Laboratory 

244 Wood Street 
Lexington, MA 02420-9108 

ABSTRACT 

The air transportation system faces a challenge in 
accommodating growing air traffic despite an inability 
to build new runways at most major airports. One 
approach to alleviating congestion is to find ways of 
using each available runway to the maximum extent 
possible without violating safety standards. Some 
decision support tools, such as the Final Approach 
Spacing Tool (FAST) that is a part of the Center 
TRACON Automation System (CTAS), are specifically 
targeted toward achieving greater runway throughput 
by reducing the average landing time interval (LTI) 
between arrivals at a given runway. In order to 
understand the potential benefits of such innovations, 
techniques for detecting spacing inefficiencies and 
estimating potential throughput improvements are 
needed. This paper demonstrates techniques for 
analyzing radar (fata from actual airport operations and 
using it to validate, calibrate, and extend analyzes of the 
FAST benefits mechanisms. The emphasis is upon 
robust statistical measures that can be produced through 
automated analysis of radar data, thus enabling large 
amounts of data to be analyzed. 

INXRQPUCTliQN 

A number of analytic and simulation studies have 
attempted to assess the potential benefits resulting from 
deployment of the Final Approach Spacing Tool 
(FAST) that is a part of the Center TRACON 
Automation System (CTAS). *• 4 One of the primary 
sources of FAST benefits is the increased precision of 
control, which is presumed to reduce the average 
landing time intervals (LTls) at each runway. In 
general, it is assumed that achieved separations contain 
some amount of excess spacing (not required by 
separation standards) and that by allowing more precise 
control, this excess is reduced in a uniform way for all 
arrivals to which FAST advisories are applied. By 
saving a few seconds of runway time for each arriving 
pair, this mechanism provides an increase in runway 


*Thii work wu performed for the NitioiuU Aeronautics and Space 
Administration under Air Force Contract No. F19628-00-C-0002. 

♦Copyright O 2001 by M.I.T. Published by the American Institute of 
Aeronautics and Astronautics, Inc., with permission. 
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capacity. The delay savings that accrue over an 
extended period of operation are found by integrating 
the delay reductions achieved over a variety of traffic 
and weather conditions. 

In this paper, data from actual airport operations is 
analyzed and applied to the problem of validating, 
calibrating, and extending the model for the key FAST 
benefits mechanism - landing time interval reduction. 
The analysis of actual operations dr* is also helpful in 
prioritizing research activities to focus upon areas 
where the greatest opportunity lies. This work extends 
the capabilities used in earlier data analysis conducted 
by Boswell and Ballin and Erzberger. *• 2 The emphasis 
is upon robust statistical measures that can be produced 
through automated analysis routines, thus enabling 
large amounts of data to be analyzed. 

PACKAGE FOR ANALYSIS OF RUNWAY 

OPERATIONS ff ABO) 

All major airports acquire and archive radar data on 
traffic in the terminal area using the Automated Radar 
Terminal System (ARTS). When combined with basic 
flight plan Information and knowledge of the runway 
layout, this data provides insight into the flow rates into 
the terminal and the manner in which particular 
runways were being utilized. A software package 
called Package for Analysis of Runway Operations 
(PARO) was written to automatically process such data 
and produce analyses relevant to the efficiency of 
operations. PARO is written in the C++ programming 
language. Data analyses presented in this paper will 
focus primarily upon analysis of four DFW data sets 
that were available during the software development 
period. These data sets were used to develop the 
analysis techniques and give some preliminary insight 
into AFAST benefits questions. Analysis of additional 
sets of data are being analyzed currently. 

PARO processing takes place in three phases 
designated GO, G1 and G2. Phase GO involves the 
reading of raw data files, correcting certain errors and 
anomalies, and producing new radar data files. In the 
original data files, tracks appear in order of the time of 
the first radar report in the track. Phase G1 involves 
reading the GO radar data files, correcting and 
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validating the input data, estimating velocities, and 
conducting certain analyses that require complete track 
data. Phase G2 involves reading and processing the 
summary data files produced from G1 processing. 
Among the variables that may be analyzed are the path 
length flown, the time of crossing the outer marker, the 
interanival interval relative to the preceding arrival, etc. 
By operating only upon the summary flies, G2 analyses 
can run more rapidly without having to process the 
more voluminous track data. 

Bavesian Rnnwav Assignment Procedure fBRAPl 

The ability to properly assign each observed operation 
to a particular runway is essential for reliable analysis 
of multi-runway operations. If radar data were 
complete and of sufficient accuracy, such assignment 
might require a simple comparison of the surface 
intercept projection of tracks with the known runway 
locations. However, several imperfections in the radar 
data (particularly altitude coverage limitations) lead to 
the need for a somewhat more sophisticated approach to 
runway assignment. 

A Bayesian approach to runway assignment has been 
developed as part of PARO. Under the Bayesian 
approach, the runway assignment is viewed both as a 
parameter that determines the likelihood of any given 
set of radar observations and as a random variable that 
has its own probability distribution. The Bayesian 
approach allows an optimum utilization of all available 
information about how runways are being used and 
what was observed with radar. The result is a runway 
assignment algorithm that is more accurate than any 


algorithm based solely upon radar data for a single 
track. 

Date ComnltteMM* 

The completeness of the data is of great concern when 
evaluating the efficiency of airport operations. Missing 
tracks create gaps in the arrival stream that can be 
mistakenly attributed to system inefficiencies. As a 
general rule, data should be approximately 99% 
complete to perform all the PARO analyses of Interest 
(That is, not more than 1 aircraft in 100 should be 
missing from the set of radar tracks). The DFW data is 
judged to be adequate in this respect 

GENERAL INSPECTION OP AIRPORT 
OPERATIONS AT DFW 

In this section we will discuss some general attributes 
of the traffic flow that are relevant to the analysis of 
imerarrival spacing. Figure 1 shows the runway layout 
at DFW. There are seven runways and thus 14 possible 
landing directions. When traffic is flowing to the north, 
the airport is said to be in a "north flow". When traffic 
is flowing to the south, the airport is said to be in a 
"south flow". 

Figure 2 shows a selection of tracks plotted during a 
period when the airport was in a north flow. It is 
difficult to determine exactly how efficiently the 
runways were being used by casual inspection of the 
actual tracks. However, the plots and analyses that will 
now be described are designed to provide insight into 
this question. 



-3 *2 *1 0 1 2 2 

x(nmi) 

Figure 1. Runway layout at Dallas/Ft Worth Internationa] Airport 
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The bxir l>FW daw .icu analyzed arc listed in Table ). 
These data sets contain over TSfX) trucks of which about 
half are arrivals. The weather for all data tettc was 
VMC with reported von In lay '5 exceeding IQ unit. 
However for data set DFW.03, a penod oflMC weather 
coded im! y lour hours before the data set began 

Figure 3 depicts the time history of operations for one 
of the DFW data sets. In till* figure, the time of each 
individual arrival and departure id tbowit In association 
with the runway of operation. Several features of the 


traffic; flow nan be seen. Nolc that «! approximately 
10:10 there is a change in runway configuration from 
”001111 How* to ‘south flow**. The rate of operations 
varies greatly with time. fNrriods of intense activity 
lasting for 45 60 minutes ire followed by lulh in which 
only modest numbers of operations txcur. Such 
irregularidcH arc attributable primarily to airline 
scheduling, but they can also be prx»duccd by the impact 
of convective weather upon traffic flow and approach 
twites. 



Table 1. DFW dalu sets analyzed. 
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Dtua fwt DFW.01 



7 e 9 10 11 12 13 

Figure $. Kunway ntiitaition showing a change in runway cotitigu ration far <kta set WW.Ol ut 1.0:30 focal 

time. 



While the pi<a.s shown ixi the previous act lion provide 
some Insight info how the airport wts operating. they 
do not atk*w us to assess with any quantiUtuvc precision 
the ctlk'fcmy of the spacing being achieved tor any 
paitkohtr runway. Itt part, tin* n; because the spacing 
achievable under radar sepantlti Hi tnandards varies with 
aircraft 'weigh*, clans, approach speed, approach 
geometry. and other lancrs, TechpkjucH for find i 
aiufclyatii wifi now be described The key feature of the 
analyst}^ in a focus upon the lauding time interval* 
i.LTU) achieved. The LTI uf. the runway is defined as 
the time separation between two ancccssivc landing. It 
i:> Uie difference between tfie time one aircraft erodes 
the run wav threshold and ihc* time the previous landing 
aitcraft erossul the same thresh'd?!. (Note: It iu 3l::o 
possible to mcaiiote 1,7 Is «t the outer marker fOMj, but 
such nme intervals will not be employed in thin paper.}. 

The ihroujrbpui of u runway over any arbitrary time 
period b fiimpiy the inverse of the average Lfl daring 
that period. lor example, if the average 1/1*1 i* 120 
seconds. then the throughput must be 1/1*0 aimT*ifi/soc 
or 30 aircraft/ bout Wc will detir-c the capamty of o 


runway ua the sustained throughput achieved under 
saturated traffic ctxtdiL ions Titen the mean LTI under 
saturated conditions is the inverse of ;hc runway 
capacity 

Figure d allows um L.TU for three hours of operaritMis at 
Dallas Ft Worth international Airport fdma set 
DFW.04). Six arrival rushes are clearly .seen during the 
10 hours of data. This plot provides 'trough*. into the 
extent to which the loading upon the arrival runways 
was balanced. It can be seen that during the $ :00AM 
push, runway 350 «*% not ns heavily ioade-J as the 
other three runways. But during the 9:30 AM push, all 
four runways appeal to have been loaded equally. 1 here 
were brief periods in which landing intervals of 60 
seconds or so were achieved for seveia! successive 
aircraft. 

W hilt w* have ct-osen to measure LTIs at Hi*? runway, 
i nferarri vai times can also tic measured at the ontcr 
marker. A comparison of the LTi measured at the outet 
marker :i c«-i the rim war fa provider! Iti Figure 5 (using 
data from data set t>FW.O \ y ft can lx- seen that ihcrc 
doe* no* appeal to be any clear tendency for the IT In 
b> become either ^rioter ;>i smaller between the outer 
marker and the runway. This implies- that control 
actionii taken within the outer marker art: not 
significatitiY chnogtn?. the tnierafnv:il limcx. 
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Time of Day (HR) 

Figure 4. Landing time intervals for dote set DFW.04 



LTI at Runway (eec) 

Figure S. Comparison of t.Ttv measured at outer niarVer and at the run wav (datii set DFW.Oln 
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A histogram o f the observed landing lime interval for 
tbo four combined data «*,u is provided in f igure 6 
The most common separation in the 90-100 second 
rang*. Figure. 7 provides o similar histogram for the 
minimum In-trail separation* observed for 801 arrival 
pairs in which both aircraft were in the “large" weight 


claw. .Separations below 2 nmi appeoi to be mostly due 
to the occasional use of visual procedures in which 
altitude separation was maintained visually. In both 
figures a line showing a theoretical ill to the hunogiain 
is shown (An explanation ot the theory follows.) 




ftitt Tim* Interval { me) 

• mkt p^.nf of UyjhiiW 'i 

W&m re. 6, Landing time interval histogram for continued UFVV arrivals (1758 pairs). 
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isa 4 -- 


1 00 



Minimum Irv Trail Separation (NMI) 

figure 7. Minimum ift4rail separation <da»crved for 801 pairs with both aircraft in tine* Urge weight claw 

(DfWj. 
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VANPEVENNE MODEL FOR INTERARRIVAL 

SEPARATIONS 

Any analysis of actual runway operations must 
recognize that at some times the traffic flow will be less 
than the runway capacity and that gaps will occur 
between aircraft that are not due to any inefficiency In 
the spacing process. A statistical model that takes this 
into account helps avoid confusing these gaps with 
excess spacing inserted by the final spacing process. 
Vandevenne 3 developed such a model for the 
distribution of observed interarrival separations. PARO 
employs the Vandevenne model to provide a more 
robust analysis of final spacing performance. This 
section describes that model. 

The Vandevenne model is motivated as follows: Let us 
assume that controllers attempt to achieve an 
interanival time separation D that represents the closest 
comfortable target spacing for specified separation 
standards and operational conditions. The actual time 
separation achieved, S, will differ from D for two 
reasons. First, there is imprecision in spacing. Second, 
there may be gaps in the arrival stream that are too 
large to be dosed by the level of control available. The 
Vandevenne model assumes that the errors and gaps are 
additive so that 


S = D + e + g (1) 

where c is the imprecision error and g is the time gap 
that cannot be closed. The model assumes that e i s 
normally distributed according to N[0, o 2 ]. 


Vandevenne noted that if the arrival stream is random 
at a given average arrival rate X, the time gaps between 
arrivals prior to application of any control actions will 
have a Poisson distribution such that 


f g (x) = Xexp(-Xx), x*0 (2) 

It should be noted that although time separations in a 
single arrival stream will not be random because of in- 
trail separation standards, the merging of multiple 
independent streams results in an initial set of 
interarrival times that is approximately Poisson. 
Vandevenne assumed that all interarrival spacings will 
include a time gap component, and that this time gap 
will have a Poisson distribution. The components of S 
are summarized in Table 2. 


Vandevenne showed that the resulting probability 
density function for S is 


f.(y) 


Xexp 




( i-p- 


2^ (3) 


where Fsh is the standard normal distribution. 


In many analyses of actual data, the value of X changes 
during the period of observation. This violates the 
assumptions in th? Vandevenne model. For that reason, 
X should not be viewed as providing a good indication 
of the actual arrival flow rate in the data. 1* is better to 
view it as merely a parameter of the distribution that is 
used to correct for the existence of time gaps in the 
interarrival time observations. 


Table 2 The Vandevenne Model 


Variable 

Definition 

Distribution 

D 

Time separation that 
controller attempts to 
achieve. 

Fixed for a given aircraft pair 

e 

Error in achieving targeted 
time separation 

normal, zero mean 

f - (x) ■ ^r xp (-5?) 

s 

Time gaps in arrival stream 
that cannot be closed by 
control in terminal area. 

Poisson 

f,(x) - Xexp(-Xx) 
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H*e targeted lime spacing. 0. is a key parameter i-iace 
tlie inverse of D Is me inherent capacity of ih^ tumvay. 
Note chat when wiftitursiHtf How exist.*; the mean 
•swerved spacing van ire significantly greater than i). 
A |XrisU ive hlas result if the mean spacing is asaumet) 
tiy be equal to the tarried spacing. The Vnndcvcmie 
nK*de! can produce a aeati y unbiased estimate of D 
under sueh conditions. This result* in a more robust 
analyse* that is better suited ui uutxmmted processing, 
foqxmenee has shown chat the form of the V and avenue 
model provides; a gcxx! fit to actual data, it ha* the 
essential characteristic: of a major peak reflecting the 
predominant normally distributed errors ard n Kmg tail 
reflecting gaps arising from other process. 

hJ! ^ 

We will now discuss how ihe ptuametets of the 
Vandcvenne distribution relate to nittway capacity and 
to potential FAST hetstfits. 

T he capacity of a runway (defined as the sustainable 
throughput when sou ik a ted with traffic) is 
approximately WD. FAST capacity benefit ace 


aMumrf ;o to derived from ivducdons in the value of 
D. 

At first glance, it appears that the parameter o tuts no 
effecr on csipadty since the spacing error it produces 
terufv to average to /-ctu, However, it is commonly 
(Guwncrfi that in actual operation* the statue of D is 
affected by o ixcause of a need to insert a safety buffer 
between eac h pair of aircraft. This buffer guarantees 
that imprecision will not cause frequent violation* of 
separation standards. The »!» «f the totter is selected 
t:> keep ihe rate: of separation violations below some 
ievd # </., ff o is decreased, the safety butter can be 
decreased. For ihe Vacidcvencte model we can model 
fhe target separation ns 

i!>«:D,, + «F 4 ..\l.<x} W; 

where t> 0 Is the required minimum lirnc separation. 
f'W ! the inverse of the standard normal diiiixilnnion 

(with acrw mean*, and a h the allowed rate of violating 
this -wparatloo. Figure fc shows how das capacity of a 
rumvs-v is affected by the value of «:i when the 
uncertainty buffer corresponds to either la 01 3a. 



O (sec) 

Figaro $. F fleet of imprecision (a) upon* un way capacity. 
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In many cases, the major determinant of D 0 is the in- 
trail wake vortex separation standard. This standard 
depends upon the aircraft weight class combination for 
a pair of successive arrivals. In translating the distance 
standard to an equivalent time standard, we must also 
consider the speed profiles of the aircraft on final 
approach. 

To provide a more relevant comparison of aircraft with 
different weight classes and speeds, we will usually 
subtract the computed separation standard from the 
observed separation to yield the excess separation S* 
defined as S* = S - D* 

When the value of S* is negative it means that the 
actual separation achieved was less than that indicated 
by the applicable radar separation standard. This does 
not necessarily mean that any standards were violated 
since under visual meteorological conditions the radar 
separation standards do not have to be applied to 
aircraft that have their traffic in sight. 


The advantage of using S* is that it allows combining 
pair separations values for all aircraft types under the 
assumption that the applicable values erf o and X are 
independent of an aircraft's weight class and final 
approach speed. This assumption appears justified by 
data analysis completed to date for Dallas/FL Worth, 
but should be verified again when different airports are 
analyzed. 


MAXIMUM LIKELIHOOD ES TIMATION OF 
MODEL PARAMETERS 


Given a set of imerarrival time separations, how do we 
go about finding the model parameters for fitting the 
Vandevenne model to the data? Vandevenne suggested 
using a maximum likelihood estimation technique. 
Suppose that we observe N arrival pairs. Let the i* pair 
have separation y,. Then the log likelihood function is 


L » 



where f y is the probability density function for the 
separation. The maximum likelihood set of parameters 
is the set that maximizes this function. Vandevenne 
found the maximum likelihood values by generating 
contours of likelihood and using search techniques on 
these contours. While this method is theoretically 
sound, the estimation of the likelihood function for each 
point on a contour involves N separate evaluations of 
the density function f r If large databases containing 
tens of thousands of arrivals are to be analyzed, the 


computational load could be a hindrance to the analysis. 
For this reason, an alternative technique was developed 
that computes an approximate likelihood value directly 
from the histogram. For this technique, the N data 
points are compiled into a histogram with H bins. The 
likelihood factor is calculated for a separation at the 
midpoint of the histogram bin. The mmg factor is then 
assumed to apply to each point in the bin. For example, 
let the count of separations falling into bin 1 be Let 

the mid-point of the separation interval for bin i be . 
Then the approximate likelihood function can be 
written 


L « ln[f y (y,)] (6) 

With this approach, the number of times the f y function 
must be computed is equal to the number of histogram 
bins instead of the number of points within those bins. 
This greatly expedites the search for the maximum 
likelihood values. Inspection c.‘ several cases indicates 
that as long as the histogram bin width is less than 
approximately one-half o , the maximum likelihood 
parameters derived in this way are almost 
indistinguishable from those derived by using all N 
original data points. 

What is the accuracy with which PARO is able to 
estimate the three parameters of the Vandevenne 
distribution? Clearly, the accuracy will depend on the 
number of points that are available for forming the 
estimate. It will also depend upon the bin size used in 
the histogram. Table 3 presents simulation results for a 
case in which the true parameter values art D = 72.0 
seconds, o- 18.0 seconds, and X » 90/HR. For each 
entry, Monte Carlo simulation was used to generate 100 
histograms, each with a bin size of 10 seconds. The 
standard deviations of the parameter estimation error 
decreases roughly as the inverse of the square root of 
the number of points used to construct the histogram - 
an expected result For o, a standard deviation of error 
that is less than 10 percent of the true value is achieved 
with 400 data points. 


Analysis of LTIs from DFW 

The LTI distributions that exist in the four DFW data 
sets were analyzed by first generating LTI histograms 
for each set separately and then for the combined data. 
The maximum likelihood fit of the Vandevenne 
distribution to each histogram was computed. Figure 9 
depicts the histogram of excess LTIs for all four data 
sets combined. 
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Table 3, Landing Time Interval tLTI) Analyst* 


r No. Pfliatt 1 
io 

iimwuin 1 

! MkwlEhm 
ia.D 

Std.Pev.of 

D 

| ldaa£tm 
feuz 

MD«r. 

uLa 

M&ro .Iks*, j 

!<L6 

I SKi.&vaf - 
1 4 - 

200 1 

•0.254 

3.249 1 

-0.4S i 

2.730 

1.954 1 

10.239^ 

400 

0.307 

2.483 

6.107 

1.78) 

2.562 1 

7.150 

800 1 

42 56 

1.774 I 

•6.216 

1.300 

0.471 | 

4.723 

1600 

-0075 1 

1.009 

4).(K)8 

0.8 tO l 

6.536 1 

3.303 


200 


ISO 


100 


50 


0 
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ah am tseu- me ?*•!*. 


Varstevenne maximum ifcoihooa fit 
D*7.n. o* 177 g f *.« 654 AC/br 


Cxceaa Lending Time Interval (see) 
fm*3*eeiit c*f 10 aftcbir^ 

figure *>. Distribution of excess kiiterarrWal spacing* (S*) for 1758 DFW arrivals 


line rtmilu of the analysis an? summarized in Table 4. 
Li can I k j««n from ihc histogram* ih at the siuipe ot the 
distribution closely resembles the Vandcvcnnc 
distribution. The maximum likelihood parametea are 
similar lor all acts. The combi oed value of p was 17,7 
seconds. Hiis U only slightly less than the 19*20 
second values reported by Ballin acni Err.bcrgex . 2 

The tan that D tends to be slightly le*s than zero means 
Uiat the strcralt were often Achieving intervals smaller 
than would be possible under radar separation 
standards. This implies for this predominantly VMO 
da*.;, it »s unlikely that. ;ia AIAST calibrated to preserve 
radar separation standards could have increased 
throughput by reducing the "imprecision buffers'* 
im>>rp*>:aied m l> Nevertheless ARST might have 


been able to provide benefits by anticipating ami 
removing the larger intcrarrival gap*: that arc related to 
the overall flow to the runways. Additional data 
analysis is being pursued to confirm this and to address 
the same question for IMC condition*. 


SE PARATIONS 

This section present the results *>f several lype* of 
analyst* conducted to investigate the reasons for the 
differences between the LT1* obtamed under different 
conditions. 
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Correlati on Analysis for Lite 

Otic way of marching for factors that aifeci final 
spacing efficiency is to examine the linear correlation 
coefficient between various variables am) the excess 
separation This ty{>c of analysis can fail fo detect 
certain types of noalincar dependencies, bat vrill 
nevertheless identify a number of relationships that 
deserve further scrutiny. 


Figure 10 shows the parameters used fio characterize the 
final approach geometry. The parameter, arc defined in 
the runway coordinate system for winch tlie origin ;s 
the runway threshold Vkr. pit ra me ter YBA8K U the y 
v^lue aL which the base leg was; established. The 
n&ntrneier lit..? ii Hie v value 2* which the ii/CfAl'l 
achieved flight ;Uwtg the cem:r\U,c <if ihc runway. The 
criicriu For “ccittedloe* status i* that the aircraft has to 
be within 2 uroi of the centecline and h&vt a heading 
within I S deeteeK of the runway Iwartmo 

^ # •• - ^ IT 


Table 4. Errors in Estimation »f Vandevenne Parameters 


Data Set 

No. of Aircraft 
Pairs 

D 

(sec) 

a 

(sec) 

>. 

(AC/br* 

DPW.01 

418 

- 7.9 

14.3 

62.8 

DFW.02 

391 

-0.3 

17.2 

69.0 

rvrsir Ait 

N 4A 


)y.v 


Ut H .U7 

4.1 * 


oe.i 

DFW.04 

707 

•12.3 

17.6 

64.4 

Ail sets 
combined 

1758 

-7.1 

U. ~ - J 

17.7 

65.4 
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Table 3 defines the “W variables’* that were used to 
describe the pair of aircraft generating a single LTI 
value. 

Table 6 provides the serial correlation coefficient p, for 
these variables when correlated against the full value of 
the LTI. Only large/large weight class pain were used 
to avoid variations due to differing wake vortex 


separation requirements. If the initial gap between 
arrivals was too large to be closed by typical control 
actions, the pair was excluded from this analysis. In this 
table, p is the probability that the observed correlation 
coefficient would be as far from zero as observed if the 
actual value were in fact zero. 


Table 5. Descriptive «W Variables” for an Arrival Pah- 


Variable. W 

Definition 

absolute LTI 

Absolute value of landing time interval (sec). 

AP_PATTERN 

Approach pattern type of follower (1000=straight-in, 2000=downwind/base) 

mmmmssi nsa 

Azimuth at which aircraft first appeared (degrees). 

LTL.OM 

Landing time interval at outer marker (sec) 

LTlmin 

Minimum LTI permitted by separation standards (sec). 

pathlength 


S* 

Excess landing time interval at runway (sec) 

t_on_CL 

Time spent “on centerline” (CL) state before landing (sec) 

V2/V1 

Speed ratio (final) of follower to lead aircraft 

vel_op2 

H <?♦ H ♦ FIB i I P» wi a j j 

vOM2 

Velocity of follower at outer marker (KT). 

WEIGHT_CL _DIF 

Weight class code of lead minus that of follower. 

wtdassjead 

Weight class code of lead (l=light, 2=large, 3=B757, 4=heavy) 

YBASE 

y coordinate of base segment for follower (meters). 

yCL2 

Centerline intercept coordinate of follower (meters) 

|HHn3E3Z35I 

Difference in centerline intercept coordinate of follower and lead (meters). 


Table 6. Linear Correlation Analysis of Foil LTI: S vs. Variable W, Large/Large Weight flags Pairs 


Variable, W 

n 

means 

Std. 

Deviation 

ofS 

■ 

Std. 

Deviation 

ofW 

P 

p 

signi- 

ficance 

V2/V1 

1211 

105.86 

42.890 

1.014 

0.148 

-0.149 

0.00000 

• •• 

vOM2 

1211 

105.86 

42.890 

86.3 

14.2 

-0.120 

0.00003 


vel op2 

1211 

105.86 

42.890 

67.221 

■H5D1 

-0.082 

0.00459 


A2 F1RST2 

1211 

105.86 

42.890 

299.6 

455.8 

-0.045 

0.11802 


vSSSSSSaKKKM 

1211 

105.86 

42.890 

24S.8 

9371.4 

0.011 



AP PATTERN i 

1211 

105.86 

42.890 

10673 

263.3 

0.020 

0.49171 


WEIGHT CL DIF 

1211 

105.86 

■KEES3 

0.015 

0.768 

■cuni] 

0.08173 

• 

wtdassl 

1211 

105.86 

42.890 

2.045 

0.609 

0.092 

0.00143 

••• 

LTlmin 

1211 

105.86 

42.890 

0.4 


hsu 

0.00001 


yCL2 

1211 

105.86 

42.890 

-17056.9 

9231.8 

Hsu 

0.00000 

••• 

pathlength2 

1211 

105.86 

■EIEJ 

150929.8 

23251.4 

0.149 

0.00000 

•IS 

YBASE2 

1043 

107.42 

44.733 

-173253 

54283 

0.157 

0.00000 

Ml 

S* 

1211 

105.86 

42.890 

29.0 

43.6 

0.905 

0.00000 

M* 


1211 

105.86 

42.890 

■EH 

45.0 

0.911 

0.00000 

• •• 


Significance code: • ■ significant at 0.10 level, •* = at 0.05 


level, ••• at 0.01 level 
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The following observations apply to Table 6: 

• Separation increases when the total path length 
flown Increases (p m 0.147). This may reflea 
greater constraints and traffic interactions 
encountered by aircraft that fly longer paths. It 
could also reflea the fact that having to maneuver 
within terminal airspace introduces imprecision 
into the spacing. 

• There is high correlation (p = 0.880) between LTI 
measured at the runway and LTI measured at outer 
marker. This suggests that if efficient spacing isn’t 
achieved at the outer marker, then it is uolikely to 
improve much due to actions taken within the 
marker. 

Table 7 provides a correlation analysis for W variables 

correlated against the excess landing time interval, S*. 

Here all weight classes can be combined. Note that 

• Excess separation is negatively correlated (p » 
-0.203) with weight class difference (lead minus 
follower). This indicates that when the lead 
aircraft is heavier, the separation relative to wake 
separation standards is less. 

• Excess separation increases when centerline 
intercept is closer to the runway (p=0.135 for 
yCL2). Excess separation decreases with more time 
spent on the centerline (p=-0.122 for L.OIUCL). 
The reason for this is not clear, but may have 


something to do with the ability to tighten 
separation through speed control as compared to 
trying to achieve tight separation by a precise turn 
from a short base leg. 

• Excess separation increases when pathlength 
increases. (See earlier comment for Table 6). 

♦ There is negative correlation (p » -0.254) with 
absolute LTI allowed by separation standards. 
This suggests that there is a tendency to space 
closer than the standard for the larger standards, 
perhaps through use of VMC procedures. 

There is high correlation (p » 0.905) with excess LTI 
measured at outer marker. Again, this indicates that 
actions taken after the outer marker have little impact 
on the final time separations. 


A 




L AL ANALYSIS OF FACTORS 
AFFECTING LTIS 


Differences Between Ronwavs 

An obvious question to ask is whether LTI distributions 
are the same for alt runways. Figure 11 provides 
histograms of LTIs for each runway at DFW using the 
combined data sets. The LTI distribution for each 
runway appears to be similar except possibly for 
runway 18L for which very few arrivals were observed. 


Tabic 7. Linear Correlation Analysis of Excess LTI: S* vs. Variable W, ail aircraft pairs 


Variable, W 

n 

mean S 

Std. 

Deviation of 
S 

mean W 

Std. 

Deviation 

ofW 

■ 

p 

signi- 

ficance 

V2/V1 

in 

101.36 

36349 

1.017 

0.154 

-0.138 

0.00011 


vel op2 

Ea 

10136 

36349 

67.672 

wmosm 

-0.123 

EESHi 


vOM2 

HFl 

10136 

36.349 

86.7 

13.6 

-0.035 

032483 


A Z FIRST2 

Em 

101.36 

36.349 

Til 

448.8 

-0.025 

0.48288 



Em 

■HF1 

36349 

167.1 

■cm 

0.035 

0.32807 


AP PATTERN 

787 

10136 

36.349 

1079.9 

281.5 

0.048 

0.17857 


LTImin 


10136 

36349 

-03 

17.6 

0.116 

0.00121 

•— 

yCL2 

787 

10136 

36349 

-16967.1 

9141.7 

0.144 

0.00006 

•M 



10136 

36.349 

1527263 

22616.6 

0.147 

0.00004 

••• 

1 YBASE2 

Ell 

102.98 

Hsni 

-17549.5 

5447.4 


0.00000 

• M 


Em 

10136 

36.349 

101.7 

37.9 

0.880 

K'l'M'MJ 

•— 

[s : 

Eli 

10136 

35.349 

31.8 

353 


BftlU'I'J 

m— 
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Figure 11. LTI Oil Terences between runways 


Rrure 12 examines the effect on the LTI distribution of 
ihe epeed (alio between (he lead and following aircraft. 
If b is (now difficult lor controllers to anucipote the 
effect of such *pc cd difference^ then differences might 
he expected to appear in tiie histogram*. For the data let 
Figure \2, the values ttf the targeted tiri^e separation. D, 
are quite similar except possibly t*x lownr D value for 
the case of a Following aircraft more itiax: 20% Faster 
than the lead aircraft In general, it appears that 
contiolicrs at DFW am quite skilled at to kit!?, the 
differing landings speeds of aircraft into account when 
spacing them. 

It seems possible that the precision of tniewrtiviJ 
sparing can be affected by the geometry available for 
mailing the final spacing adjustment. For aircraft shut 
fly n downwind segment, (he controller i;» able to 
choose the hxati or of live bm,* segment u* achieve 
proper spacing. Hu? ‘he mm required doing so van be a 
source of imprecision. Ik there a difference in spacing 


performance that should be addressed by foots »ncb as 
HAST? 

An analysis of approach patterns wo* conducted fi> 
dividing all pair* of successive aircraft into lour groups 
depending u|x»n whether or not the approach Involves! a 
downwind base trajectory. Approaches without a 
downwind phase were called "straighten*', aUlumgh it 
should be noted that *xne of ihe.se aircraft approached 
the runway centerline a! a fairly large angle. The sei 
labeled ''downwind >3uai?hr includes nil pairs for 
which an aironift on straight -io approach war. followed 
by an aircraft flying a downwind segment The font 
histograms that result arc shown in Ftrurc 13. A( 
DFW, aircraft are. generally first directed io the 
coffierposl fixe* that are nKtft consistent with the 
north/sovth direction of flow, and hence straight-in 
paittms predominate for final approach. While there 
arc no statistical!? significant differences in the S f 
diaUxbulion for this data, it d«*r> appear that aircraft 
following a downwind feodci lend tv have higher 
media a values 
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One aiidiiicnnl question concerns the relationship 
between time intervals amt actual In trail spacing. i» 
this re$?ott. we have expressed scporulicms in terms of 
time, but actual separation standard* are expressed in 
distance. Figure 14 show? the rctelionslup between 
excess landing time Interval* and excess in-trad 
spacing* for arrival pairs in data set DFW.GI. It can tx 
seen that there is high correlation cp « 0 X43j between 


the two xncanetc*. The doiied line shows the space to 
time conversion that would if tin? excess spacing is 
traversed by the trailing, aircraft ul a speed of L50 kiKHS 
This analysis suggests that conclusions about 
performance deduced from Inspection of time 
separation* are likely to be Uk same tx* those of an 
analysis that used uniy ejXUiuI vcptm lions. 


tot OFWjOi 



Excess Interarrival OUiunce (NMi) 

Figure id. Comparison of excecv spacing expressed in tiou< and distance 

Nor is it clear that they can be sustained in periods 
of prolonged saturation. 

To general, the LTI achieved at the outer marker is 
preserved at the runway. Variations appear to be 
random with no discernible tendency to change in * 
given manner. Thus, them is httie evident that 
visual separation practices applied within the outer 
marker are having a sigitific&ni Impact upon 
spacing?. 

1 he targeted separation. appears to be uixiui ?2 
seconds, which corresponds to a throughput of 
about 50 aircraft/hour. The fact that this rate is 
almost never austained in practice suggcsti; tha» 
there is art opportunity lo increase dmxidiput if the 
consistency oi flow to the final vector position U 
improved. 


VmS&MLQZ&J 

This work has developed robust, computationally 
practical data analyuts routines Uu*l can lx used to 
provide insights into runway operations through 
analyst:-: of radur data. The questions that can be 
addressed am, relevant to the benefits mechanism* of 
Af-AbT and to determining the total benefits that can Ik 
achieved with A FA ST imptemenuuion. 

11m following observation*: apply to Uw four data sets 
analyzed for DHV It should be carted that all these 
4iU» sets involves* VMC weather and the generniKy of 
the cunriusuaw drawn fcom this- limited data set ha*, wm 
txxt: proven. 


In short rush periods {of Ml minutes or &>;, very 
high landing rates of 60 AC/Hfc or more tire 
obtained on single runways It in not clear that the 
high pcai: rates observed enn lx* obtained in CMC. 
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• Targeted separation (D) tends to be slightly less 
(by about 7 seconds or 0.3 nmi) than the value that 
would be expected if radar separation standards 
were the sole determinant of target separation. This 
suggests that visual procedures near the runway 
may have allowed separations to be tightened 
enough to overcome the effects of any “safety 
buffers’* that were applied during the earlier radar 
separation process. 

• The occasional presence of larger LTIs during 
periods of saturated flow is a further Indication that 
irregularities in the flow may be contributing to a 
loss of throughput in VMC. This phenomenon 
deserves further study since the ability of AFAST 
to reduce such irregularities provide capacity 
benefits under VMC conditions. 
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ASSESSING DELAY BENEFITS OF THE FINAL APPROACH SPACING TOOL (FAST)* 

terry D, WuJvb and John W, Andrews John E. R->bioson MI 

MU' Lincoln Laboratory NASA Arms tfaswrr.h Center 

244 Woo j Street Htffitt F<eld, CA $403* -MOO 


htxtognm. Ml 02*20-9 JOS 
ABSTRACT 

Air itaMIc: delays aye costly, NASA is developing the 
Final Approach Spacing Tool (F AST) to help increase 
throughput and reduce tho approach component of 
nirluxTio delay, Analysis and field trials have suggested 
lhal FAST can help controllers increase arrival 
throughput on busy runways by several aircraft per 
hour. 

Two major simulation studies Imve predicted thai tteiay 
reduction? ft cm such throughput increases would «ave 
several hundred million doilara annually. The studies 
mo provided useful date on operations at major 
airports. However, ‘heir predictions disagree on delay 
savings tor some airports and omit other airports of 
interest. Their predicted delay savings For some 
airports are higher tttnn actual reported delays tor those 
airports. Because of resource and rime limitations 
neither study considered stotm disruptions to arrival 
routes, and neither addressed downstream delay 
propagation caused by schedule disruption. Both of 
these effects change the dollar savings from FAST. 
Ah hough delay propagation can only multiply delay 
savings, the effect of a storm can be positive or 
negative. Route disruption* from storms can 
temporarily pre-vein FAST from operating. But storms 
can abc leave large queues that magnify the value of 
incremental capacity increase;? from FAST when it 
returns i<> operation. 

In this paper we summarize and compare the two 
benefft models We present a simple benefit model tor 
lodging the accuracy of the two models and tor Helping 
to ranking benefits among airports. We examine 
meaiKKcd delay data published alter die completion of 
the studies that helps to validate the model results. We 
use this data to examine delay coupling between 
airports and we use it with a published model for delay 
propagation to estimate downstream delay 
nmlOplioiiitou for Dallas fori Worth Airport. Finally* 
we examine bow storms affect delay and how they 
might modify ihc estimated FAST benefits. 

*ih« wor* was fwifcrir.cd for th-- National Ace&ruiutfcs and Space 
Adrr.iiRvff alien under Air fmcc CX? :5, nictNe. £J’#2$*£0-C*0002. 

tC:>pyng?:l v 200 1 fcy M i.l* by the, Awene&n If?Miiv*ue of 

A*!CH3VU£$ ^Tid /vjl'OUSUi&S, Uu.., w;«Ji pCKtuVijcn 


The cost of aft traffic delay grows each year a« (Mavis 
and demand Inci’case and as fuel costs incteasc. The 
Air Transport Association ciifim&tcd that airborne 
delays cost US air carriers about S800M in 199*).* 

Figure l illustrates airborne delay during congested 
periods at DFW on Monday 10 April 2000. Hie figure 
shows the arrival rate at ; S-mi flute intervals and the 
corresponding average airborne delay tor arrivals in 
each 1 5 -minute period. Toe arrival rate was determined 
from radar data, and the average airborne delay was 
obtained from the f AA's Consolidated Operations and 
Delay Analysis System (CODAS).* Six distinct arrival 
rushes are evident, and associated with each rush its a 
transient build-up of queuing delay. 



UtKO (UU-ftf >.<sw) 


Figure 1 Arrival Kate and Airhorxte Delay at Dallas 
on 10 April 2000. 

As part of the Center- TKACON Automation System 
(CTAS), NASA h developing the Final Approach 
Spacing Too) (FAST) to help reduce this kind of delay .* 
FAST is intended to increase terminal throughput and 
reduce the approach component of airborne delay- 
caused by queuing congestion at the airport, U 
accomplishes this by providing planning advisories for 
efficient runway balancing and arrival sequencing; and 
by helping to increa-ie the accuracy of final approach 
spacing. 

Analyses, vimuluiions, and field trials inchoate that 
FAST could help controller* increase arrival throughput 
on a busy runway by wvetui aircraft per hour. Major 
independent studies by Seagull Technology. Inc. and 
Logistics Management Institute (1..M.1; have estimated 
the potential dollar saving* from such throughput 
improvement:*? su 10 nu\ior US airports. 4 '*** ‘ ,,v Both 
studies determined from initial analysis that FAS) has 
the potential to decrease auc.taft i? -.ter -arrival tin:^s by 
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about \5 seconds and thereby increase throughput by 
about 4 arrival* per hour per runway. This h 
accomplished by helping controller!! reduce the 
variance of inter-arrival timing. Hie two studies 
assumed thru tbia more precise arrival timing allows 
controllers to reduce imefxtrrival spacing without 
increasing the incidence of separation viol alums. They 
then analyzed the extern to which such a throughput 
increase would reduce approach delay from congestion 
during arrival rushes. 

The studies employed independent demand and 
capacity estimates and used separate queuing engines to 
calculate the reduction in queuing delay, the modeling 
process is illustrated in Figure 2. A queueing engine 
converts time-varying airport capacity and demand to 
delay. Seagull Technology used an FAA-provided 
discrete event simulation to convert capacity and 
demand profiles to delay, LMl integrated the 
Kolmogorov queuing equations to derive delay piufilea. 
They used independent, but similar, models of direct 
operating cost to aggregate and convert aircraft arrival 
delay savings to dollars. 


nwaswemam.* for the airports studied. Their models 
iiKimt) on delays In. IMC and VMC and did net reflect 
tbc fact that FAST cannot predict flight trajectories and 
mur>t cease operating when hazardous weather disrupts 
arrival routes. It fellows that they also did not account 
for the large residual queue*; that must be serviced when 
routes re-open after storms. 

The dollar savings cfUtniafes of both were kept 
Intentionally coneer/aliYe by focusing only on direct 
operating cost xaviuga. Downstream delay propagation 
costs* which can be large for days with unusually 
adverse weather nr high demand, were omitted for lack 
of a cost m«x!cl applicable to all airports. 

The ton airporo modeled by LMl were all included in 
the Seagull Technology study. Tlie overall cost savings 
predicted for Ihc ten airports by Seagull Technology 
v&s $373 M/year. LMl estimated a savings of 
ST 1 9M/yeer. The dollar estimates for some airports 
differ more signHteantly than do these overall numbers. 
Figure 3 compares the results tor tho 10 airports that 
were studied in common. 
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Figure 2. Comparison of the Seagull Technology' 
(SGT) acid Logistics Management Institute (LMl) 
Benelit Models. 

The models rut* capaeity/demand scenarios for each 
airport with today's runway capacity* and then repeat 
each scenario with the increased capacity ihm would 
result from reduced inter- arrival rimes achieved with 
the help of FAST. Seagull Technology modeled only ?. 
capacious for cacti airport, whereas the LMl capacity 
model accounted fot all major runway configuration* at 
each import. The Seagull Technology demand model 
used a single demand prctfllo for each airport, whereas 
the LMl model used 6 demand profiles per airport to 
account for seasonal differences and differences 
between sveefcday and weekend schedules. The Seagull 
Technology study examined 35 airports. The LMl 
study analyze 10 airports. 

Both studies adopted similar limitations of scope. They 
did uc‘j validate their predictions with actual delay 
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Figure 3, Comparison of Seagull Technology and 
logistics Management Institute Benefit F^sti mates. 


Ideally, the points for all of die airports would fall near 
the line of unit slope included in the figure. Hie 
disagreement, in large for Chicago* and Atlanta. It is 
very large for LaOnatdia and Dallas/Fort Worth. These 
difference* make it difficult to determine relative 
benefit nmfcmgff for the purports. 


ASasmjR luis jg sG JBl 1LS 

To help resolve these difference:; wo developed a 
simple benefit model based <m steady state queuing. 
These benefit analyses reduce ultimately to calculations 
of the change in average annual delay that results Iran 
A given capacity increase at all runways. Although 
steady state queueing lhtx>cv does not allow esuct delay 
calculations when capacity and demand vary with time, 
it does provide insight ime average annual delay trend*. 

H we assume Utah on average, atl airport runways are 
loaded with the 3cune ratio of demand so capacity, and if 
the capacity of all runways is increased identically, then 
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to first order, steady -state queuing tells us that the delay 
reduction at each airport wilt be proportional to ls :i /R. 
where N h the traffic count aiKt R it* the total number of 
iimways a\ the airport. 


This prediction was tested against the results of bofit 
randies. Figure 4 shows site LMi annual savings 
estimate for each of the 10 alr|>ort$ plotted as a function 
of NtK. Ttie data label for each airport includes in 
parentheses the number of runways capable of handling 
commercial flights ar ihar aitpott. 

The LNIl benefit numbers follow the N ? /R trend 
reasonably well. Also plotted is a linutir least-squares 
ill \o. the savings tor the Id airports* which can he used 
to estimate savings fox other Airports based on their 
operation* and runway counts. We estimated the 
savings For the FAA Free Flight Phase* I. airport* that 
were not included in the LMf. study (Philadelphia, 
Charlotte, Denver, Miami Minneapolis/St. Paul and St. 
Louie). If the LM.i savings estimate* for FAST at the 
10 study airport* ate correct, the benefits for the six 
additional airports range from SUM pst year fox 
Denver \o S'iSMl pet year for Miami. 
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figure 4, LMl FASH savings «»rtxxiate$ vs. LN Vft), 


Overall, the result.* Indicate that airlines using the lour 
largest airport* (LAX. AIL. DFW. and ORD) would 
save the most Front FAST, IF controllers using FAST 
could it; deed increase tunway arrival throughput by -1 
aircraft p<a hour at all 10 airports, airline operators 
would recover an estimated &460M annually in direct 
opvraupa costs. 


fciu 

Wv now examine the use of recently available delay 
measure*) tern* to validate the model results;. Wltcti the 
two benefits studies began, there w;« no available delay 
metric commensurate with the queuing delay 
considered in the studies. For a number of years the 
mam quaromtive delay reports have been those based 
on FA a OPSNFT data, which counts High is with 
schedule delays exceeding 15 minutes. In 1977 the 
FA A b egar to make available statistics on delays of all 


magnitude* os part of the Consoltd;tted Operations axtd 
Delay Analysis System (CODASV 


CODAS delay is gathered and archived for more than 
100 airport* and provided on the internet for authorized 
user*. The CODAS database includes delays from 
several phases of flight. Its “arrival delay and 
“airborne** delay estimates appear to be ntou relevant to 
the estimation of FAST benefit*;. Both types of delay 
are averaged every fit been minutes and reported in units 
of minutes per arrival. ft*. deriving these averages, the 
CODAS preceding algorithms count early arrivals as 
having aero delay father than negative delay. Figure S 
illustrates the differences between these delay metrics. 


Conacrffdated Operations and Delay Analyst* System 
(CODAS) 





• k*io<* I 


: * : Arttti- 


Osf «iM*9t : 


Figure 5 Definition of CODAS Airborne and 
Arrival delays. 


Airborne delay is measured relative to the Flight 
duration predicted at the time of departure. It is the 
actual flight duration (wbx&ls-up to wheela-dowu) 
minus the predicted flight duration. Airborne delay 
docs not include departure delays. The direct opera tog 
cost «>f ait borne delay can be readily calculated. 
Although some airborne delay can be caused by err 
route weather arid traffic ft ov problem*. normally one 
of its largest components is the terminal queuing delay 
that runway capacity improvements from FAST are 
intended to reduce 


Arrival delay la measured relative to scheduled arrival 
time. It is tire actual gate arrival time minus the final 
airline Computer Reset various System scheduled gate 
arrival lime for the Right. If the High-: duration 
predicted on take-oiT h the some a* the scheduled flight 
duration, the arrival delay U da? sum of the departure 
delay, the airborne delay, and the taxi-in delay. Arrival 
delay is. usually larger than airborne delay. 

At rival delay is not related to operating cost in a simple 
way. It includes the delay that results when the traffic 
flow management process hoick aircraft ox* iho ground 
to minimize airborne delay. But reducing airborne 
delay by holding aircraft on the ground does not change 
the fact that aircraft laud behind schedule. These 
schedule delays can result in downstream ripple co-ms 
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that arc more difficult to account for than the direct 
operating coils associated with airborne delay*. An 
cxaromalion of arrival delay can help determine the frrii 
benefit to be expected from an improvement in airport 
capacity. 


MQQELJ8.ESULTS AM) CODAS DATA 


Figure 6 compares the 1997 CODAS average annua! 
airborne delay at ton airports with die LMi model 
estimate* for annual delay that: would have been saved 
in 1997 by FAST and IMA at those airports. The 
ilgure shows that the general trends for the LMI savings 
estimates and die reported delays are similar. However, 
the LMI estimates for three of the airports would appear 
to be iliogically latge in that the delay savings estimates 
from a small incremental improvement in runway 
capacity exceed the total reported airborne delay at 
those ah ports. 
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Figure ii CODAS airborne delay in 1997 vs. LMI 
estimate of delay savings from FAST*. 


The Seagull Technology study showed a s imiim trend. 
It estimated that five airports would experience delay 
reductions from FAST that exceeded their CODAS 
airborne delays, with savings estimates for LAX, i/JA, 
and ORD exceeding CODAS delays by huge factors. 


Systematic underestimation of airborne delay by 
CODAS can make delay savings estimates appear high 
for all airports. Airborne delays art? computed relative 
to operator estimates of Might time ;U inkcoff. When 
aircraft operator* predict flight durations based on mean 
historical flight times; rattier than the shortest feasible 
Might times. CODAS underestimates atrlwrae delay. 

This discrepancy is moat significant fix Los Angelo. 
Roth the Seaguli T echnology and LMI studies predicted 
a saving-; of over 2.5 million rninnott/year at LAX, 
whereas the reported airborne delay at LAX was law 
titan I million minutes per year. This Is likely caused 
by traffic flow management procedures. When the- nit 
traffic management system uses ground holds and taxi- 
out delays to regulate anti meter I ho flow of aircraft to 
an airport, airborne delay Is lower them the delay 


calculated by queuing models assuming random 
activate, &mi there is a net increase in arrival delay. 

CODAS statistics for LAX coruiteumuty allow a 
restively high ratio of arrival delay relative to airborne 
delay. In 1997, the average arrival delay al LAX was 
5.9 times larger than the average airbtime dday. The 
average arrival/airborne delay ratio for the 30 busiest 
airports was 3.4, and only two of those 30 airports 
f Phoenix and Washington Dulles) had higher ratios 
than LAX. At the other extreme, LG A and EWK, 
which He farthest above the line of unit .slope in Figure 
6 f had lower than average ratios' 2.1 and 2.4 
respectively. 

Ft appears that airborne delay at LAX may be less than 
at other major airports because rite air traffic 
management system is both, motivated and able to 
regulate and meter ext route How to LAX with unusual 
consistency, IAX k unique among major US airports 
in the complexity of its airspace and the uniformity of 
its weather. Its dense traffic and conslmnod app roach 
routes provide the motivation for consistent metering to 
avoid the need for airborne holding. Its freedom from 
storms makes it possible to meter with eomUlency, 
Smoothing of the c« route arrival flow further reduces* 
unpredictability from queuing contention at runways. 
Schedules for ilighta to LAX can absorb repeatable 
delays caused by metering and thereby trade longer 
average flight times for increased schedule 
predictability. 

The high schedule predictability at LAX k evident in 
the CGi)A$ statistics- Table 1 aurmnariees arrival 
delay statistics for seven airports in 1.977, 


Table L CODAS annual arrival delay 
(minutci/ar rival) for seven airports tn 1997, 



Alt. 

nos 

?mv 

*rw< 

MX 

1XVA 

mi. 

Me*n 

14.5 

U.9 

10.9 

i </> 

in 

m 

H.0 

Std Iltv 

16 

uy 

19 

l?.:i 

f* 7 

rx 

6.5 


In 1997 the standard deviation of CODAS arrival delay 
at LAX was only 51% of the mean. The percentages 
for DOS, DFW, LV/R, and LG A wore respectively 
86%, 79%, ?8%. and 72%. Increasing schedule 
ptcdiciahlity at LAX teduoes airline cost by decreasing 
the propagation of schedule delay a to downstream 
airports. Phis allows more efficient use of crews, 
ground facilities, and aircraft. We discuss delay 
propagation more frilly below 

CODAS delay statistics provide vauable insight into 
operational differences among airports. They are also 
useful for nxM&i.itng changes m delay, correlating delay 
between airjiorts, and studying \ht effects of weather or: 
delay, os shovvo in the following sections. 
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FAST is intended to help improve airport capacity. 
Transitions from VMC Jo IMC cause mcaswable 
statistical changes m airport capacity, Therefore, 
quantifying the relationship between local 
meteorological conditions and measured delay (Le., 
using i he transition from IMC to VMC as m analytical 
surrogate for a capacity increase) can provide baseline 
comparisons for capacity modeling results. 

The PAA’k CODAS delay database includes local 
ceiling, visibility, and wind as well as* a meteorological 
condition indicator that switches from IMC to VMC 
when visual approaches are allowed at each airport. 
We used this database to examine the dependence of 
actual delay data on local meteorological conditions at 
key airports, 

CODAS defines Visual Meteorological Conditions as 
the combination of ceiling and visibility for which 
visual approaches arc allowed. To support visual 
approaches the ceiling must be 500 ft above the 
minimum vectoring altitude, which is determined by 
airport elevation* terrain clearance, and other locai 
factors. Thus CODAS VMC corresponds to “high 
Visual Flight Rules (VFR)*\ At Boston CBOS), visual 
approaches arc permitted when the culling exceeds 
2500ft and the visibility exceeds 5mi. At Dallas Fort 
Worth Airport (DFW). visual approaches ate permitted 
only for ceiling*, above 3500ft and visibility greater 
than Snii. Lower ceiling and visibility conditions are 
considered IMC. l*hat is, CODAS IMC corresponds to 
*k>w VFR*' and below, CODAS weather data come in 
either 15-minuie or hourly summaries, Any hour with 
one or more 13* minute intervals of IMC b considered 
to be an IMF' hour. In our analysis, any d«iy with one ox 
more IMC: hours between 6AM an* midnight b 
considered to be an IMC day. 

Using these definitions we find that at OF W in calendar 
year 1997, 35% of the days had one or more IMC hours 
between 6 am and midnight In act the top >3 delay 
day.*:* wetc all IMC da>v» and <8 of the top 40 delay days 
were IMC d ays. Figure 7 compares CODAS airborne 
delay on IMC and VMC' days at DFW in 1997. 

On VMC days, the mean wtui 19 minutes of delay per 
aircraft: the standard deviation wo*: 0.83 minutes per 
aircraft, and the delay on the worst VMC day averaged 
5.o mmuies of delay per aircraft. On IMC days, the 
imauiA, standard deviations, and peak delays were 2 to 3 
times larger than on VMC days. 7 he CODAS arrival 
delay at DFW in 1997 was 1 to 4 times larger than Iho 
airtioroe delay by ail atari;;!!:;*! measures, and showed a 
similar factor of 2-3 increase in IMC. 
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Figure 7. Distribution of CODAS airborne delay on 
VMC ami IMC days * DFW 1997. 

Figure $ separates the CODAS airborne delay data into 
IMC and VMC' component* for 7 important airports 
(Atlanta (ATL), DFW, LGA, BOS, Philadelphia (FHL). 
Newark (KWR), mid LAX) with varying operational 
characteristic*. 
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Figure H, CODAS airborne delay far seven airports 
oa VMC and IMC days - 1*197, 

Results similar to DFW were found for all of these 
airports: on IMC days the mean*, standard deviations, 
and peak delay* were significantly larger than their 
valuer* on VMC day*. The observation that BOS mid 
DFW delays were equally sensitive to IMC it 
wnnewhat unexpected. The sensitivity of arrival 
runway capacity to meteorologic.il condition* differ© 
significantly bttwexsn these two airports. In vmc 
reconfiguration situation*. Boston b arrival runway 
capacity can drop by nearly <0% in IMC, whereas the 
bizgest IMC arrival runway capacity drop possible at 
DFW in 197? wo* about 33%. 

Although, cm average, DF W has excess capacity that is 
not strongly influenced by reduced toiling aiul 
visibility, DFW experiences hubbtng peek* oat h day 
that temporarily approach the available VMC runway 
capacity. During these rushes, a small decrease In 
either en route i* terminal capacity can cause a Urge 
increase in delay. In VMC the queues that build up in 
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these brief period* of excess demand are quickly 
cleared after the demand subsides h tstkess lunger to 
cleat these queues in IMC, In 1997 before the new 
DFW runway became operational, arrival capacity 
could often be reduced by loss of a diagonal runway, 
resulting in even larger queues during transient arrival 
rushes and longer residual recovery periods after the 
rushes subsided. 

The tendency for airports ro experience larger delay 
means and variances on IMG days than on VMC days 
seems ro support rhe notion of local causality. Thar is, 
if we can increase LMC arrival capacity at an airport, 
we should also reduce delay* at that airport However, 
when we examine the correlation between delays ai 
airport pairs we find evidence of systematic effects 
correlating delays over the region for CODAS airborne 
delays as well a* arrival delays. This occurs even on 
mixed days when one airport experiences some IMC 
and the other exponent** solid VMC. We also see that 
correlation decreases as the geographical separation 
between airports increases. 

Figure 9 shows die correlation between CODAS arrival 
delays at EWR and L6A for all days in 1997 for all 
four combinations of meteorological conditions. The 
correlation between CODAS arrival delays on the days 
in which l.MC prevailed at both airports was 0.83. 
Correlation was equally strong on those few “mixed” 
days when it was IMC «t one airport and VMC at the 
other. The weakest correlation was for the majority of 
days when the weather was clear at both airports. Even 
on these days the correlation was significant and 
positive at 0.49. 
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Figure 9. EWR and LG A - daily CODAS arrival 
deluy corrriilioo - 1997 all day*. 

Strong correlation between EWR and LGA delays 
might be expected because their traffic is managed hy a 
common IRACON, the airport* arc close fo each other 
geographically, and they share common arrival and 



departure fixes.". Because of dm physical proximity, 
weather conditions arc ako correlated between the two 
airports; in 1997 there were only .14 days - split 15/19- 
in which one airport experienced aome IMC end the 
other experienced solid VMC; there were 1.7! day* 
when both experienced FMC; and there were 208 days 
when both experienced $<4hi VMC. 


Strong correlation also occurs between CODAS arrival 
delays at other airport pairs since delay* relative to 
schedule are influenced by the connectivity of the air 
transport network, Whitt is ^uprising Is that strong 
positive correlation also occurs between airborne delays 
ai some airport: pair*. 


Figure 10 shows the correlation between EWR and 
PHL for CODAS airborne delay. There was significant 
correlation even though these were not schedule-related 
delays, the distance between the airports is greater, and 
the terminal air traffic is managed by different facilities. 
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Figure 10. EWR and PHL - daily CODAS airborne 
delay correlation - 1997 ait days. 


Pigurc 1J summarizes the correlation coefficients 
between selected airport pairs for CODAS airborne and 
arrival delays for all dayo »n 1997. The correlation 
generally decreased as the distance between the airport 
palm increased. 
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Figure 11. Correlation coefficients for CODAS delay 
for nine airport pairs - 1997 all days. 
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Hie annual airborne delays were oof correlated for 
widely separated airports, Howcvei, the annual arrival 
delays showed small positive correlation between all of 
these airport pairs because of downstream schedule 
impacts on high delay days. For example, schedule- 
based arrival delays at DFW and ORD were positively 
correlated, probably because both airports are major 
hubs for American Airlines. 

Downstream delay caused by schedule connectivity can 
multiply the cost of large delay events. Late arrivals 
propagate through airline schedules and result to 
additional downstream dciays. Inis effect multiplies 
the dollar benefits from reduaions in initial delay. We 
estimated the magnitude of the downstream arrival 
delay a i DFW in 1997 bused on a 1998 analysis by 
Beatty c.i al of downstream schedule delay resulting 
from 500 delayed departures from DFW. 12 That study 
showed that the number of minutes of downstream 
delay resulting from each initial delayed departure in 
related to the magnitude and the departure time of the 
initial delayed flight The reUuinnuhip was modeled in 
the form DM m l+S*F>D< where DM h the delay 
multiplier, i>D is me magnitude of the initial departure 
<Ma y, and the departure' time factor, S, wan obtained by 
a linear least-squares fit to the delay data, S h greatesf 
when the initially delayed flight takes off early in the 
day. The fitted data in Beatty** final delay multiplier 
I able shows that, to first »>f<ler % S decreases linearly as 
the departure time of the initial delayed flight increases 
from 6 to 22 hours, and tirut S is independent of the 
magnitude of the initial departure delay. Figure U 
plots S versus the departure time of the initial delayed 
flight for three values of initial delay. 
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Figure 12. Departure-time factor S as u function of 
the lone of the initial departure delay. 

One can use the dday mull ip her mode) io obtain a 
rough estimate of downstream delay tor DFW in 1997 
based on CODAS data Because the CODAS database 
focuses on delay for flights to each airport and does not 
provide statistics on departure delay from an airport. 
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one must us* the average daily CODAS arrival delay 
into DPW to approximate the a vet age dally departure 
delay out of Z)FW . Tnis approximation ignores the lint 
level of delay branching that occurs at DFW itself. It 
thus umlercvtiniatec downstream delays resulting from 
large arrival delays and provides a lower bound on 
downstream delay. 

We estimated the delay multiplier for each day by 
further assuming uniform delay between 6 and 22 hours 
equal to the daily average for all arrivals. The linearity 
of S then allows us to use the tne«m value of S over that 
period (which is 0. 77/hour} io calculate downstream 
delay 

At DFW in 1997 the daily CODAS arrival delay 
averages ranged from a low of about 2 minute* per 
flight to a high of about 58 minutes pet flight. The 
uniform delay approximation is obviously well justified 
for days with very low average delay. Thu uniform 
delay approxbrnifion also appears to be a reasonable 
assumption for days with high average delay, although 
k cat) underestimate the downstream multiplier when a 
short period of very high delay occurs early in a day. 
Hourly arrival delay’s of $ or 4 hours occurred 
occasionally at DFW in 1997. However, the hourly rate 
only exceeded 6(1 minutes about one percent of the 
time, Most days with average daily delays exceeding 
30 minutes involve episodes of relatively constant high 
hourly delay distributed dwou^heut the day. 

fhu uniform delay approximation is reasonably well 
justified when calculating the annual downstream delay 
contribution of more normal delay days. At DFW in 
1997 the annual CODAS arrival delay average was 
about 10 minutes per flight and the annual standard 
deviation was about $ minutes pet flight. Although 
one-sigma days can experience brief periods of high 
delay, the occurrence times of those delays tend to he 
uniformly distributed when considering the Urge 
number of such days in a yea* 

The*# approximations applied to COD AN arrival delay 
provide a lower bound on the annual downstream delay 
multiplier of about J 2 ai.DF\Y in 1997, The stacked 
column chart of Figure 13 shows the initial delay and 
ihc lower bound on the downstream delay estimates 
accumulated for all days In 1997. 

Ilur chart /lather breaks down delays for the 94 days in 
.1997 j/iat experienced thunderstorms within 5b nautical 
miles of DFW. The total cumulative CODAS arrival 
delay for flights into DFW in 1977 was M3 million 
minutes on those thunderstorm day* and 229 million 
minutes on the remaining non -thunderstorm days. 
Because etorm days iuni larger initial delay, they also 
had larger dewnsiresm delay. Thus, the effective delay 


25 


ALAA Guttantt, & Control Conference* Montreal* Quebec* Au^u*t 2001 


multiplier was about 13 fat storm day* compered lo 
1.13 tor non-storm day*. 



figure 13* Initial and downstcam arrival delay for 
days with and without storms im DFW in 1907. 


F AST is currently unable to predict flight trajectories 
when storms disrupt arrival routes. Thus, 
thunderstorms reduce the amount of time that F AST 
can be used. However* such route disruption* ate 
infrequent* ami the benefit of extra runway capacity 
increases dispropon ioruttt* ly when the storm has passed 
and controllers must clear out residual storm queue*. 

To determine the net effect of tJumdctstonrns on FAST 
benefits it is necessary to quantify the relationship 
between hazardous weather and delay. We examined 
hazardous wcatixt delays at DFW in 1997 and at £W'R 
in 1999.’* Weekly report logs front the integrated 
Terminal Weather System (1TWS) at DFW indicate that 
there were 94 day* that bad thunderstorms within 50 
nautical miles of DFW. 14 The DFW T&ACONf logs 
chow titat on about 50 of these days the storms involved 
enough disruption to air traffic to cause delays. At 
EWK there were 34 days with thunderstorms within 
100 NM of the aiipott that caused major delays. These 
numbers are higher than the number of day n in which 
thunder storms wetc officially reported ai DFW and 
EWR. lower personnel report thunderstorms at an 
airport when they detect lightning or thunder. On 
average that occurs 4:1 days a year cl DFW and 26 days 
a year at EWR.. 

Figure 14 is a plot of the CODAS airborne delay at 
DFW on the 50 worst delay days in 1997 sorted by 
delay magnitude. The M worst days all haul 
thunderstorm activity Thirty-four of the 40 worst 
airborne delay days were thunderstorm days. Large 
airborne delays are strongly associated with 
thunderstorms. Yet* because there were many more 
storm-free days In the year, the total annual delay on 
tHorn- free days was afx.»ut 42% larger. 

Figure 15 shows the cumulative 1997 CODAS airborne 
delay separately for thundemorm d;iy. : i and ait other 


days at DFW sorted in order of descending airborne 
delay. We multiplied the average delay on cadi day by 
that day'* arrival count to otrtain the cumulative aircraft 
delay minutes. The cumulative annual CODAS 
airborne delay on thunderstorm days wa* about 415*090 
minutes. The cumulative delay on non-storm days was 
591,000 minute*. 



OJ*? J&ak 

Figure 14 CODAS airborne delay on 50 worst days 
-DFW 1997* V 

The direct operating cost to airlines at DFW in 199? 
can be estimated by multiplying tlx airlxime Delay* by 
the $19/mimite estimate obtained from the Seagull and 
LMi benefit analyses. The results total S1L2M for 
mm- storm day* and $7.9M for thunderstorm days. 



fcy V>*ity A ttolvf 

Figure 15. Cumulative CODAS airborne delay on 
days with ar.d without storrov- DFW' 1997. 

The larger coat for nun-storm days occurred partly 
because then.’ wore more day's and, to u lesser extent* 
hecau.se there warn more arrivals *>u those days. The 
cumulative minutes and dollars for thunderstorm days 
would he larger if site calculation included nominal 
delay and dollar equivalents for each cancelled (light. 
As shown above, a complete cost accounting for 
downstream delay would also increase storm-related 
costs relative to costs on non-storm days because larger 
delays cause larger downstream rippte effects 

Thunderstorms and IMC were the main contributor* to 
large CODAS airborne delays at DFW. But high wind* 
ahme were found lo cause significant delay at EWR. 11 
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Although the predominantly North-South orientation of 
the DFW runways makes it potentially vulnerable to 
crosswinds, DFW had only one day in 1997 that was 
free of iMC and thunderstorm activity but that had 
CODAS airborne delays greater than the average for an 
FMC day Cm December 9, the delay built up during 
five boar* of 20- to 2S-kt cronwinds after 1PM, but a 
long frerxod of delay in the morning when tltc winds 
were below 10 kt aiso contributed to the high daily 
delay. Unlike EWR in 199$, where winds alone caused 
numerous large delay events, DFW in 1997 did not 
experience significant delay contributions from high 
winds. 

Delays can also be caused by inefficient handling of 
arrival traffic or by contention for air space and 
runways in peak arrival periods. But, we found that 
days with high average delay at DFW' have statistically 
lower daily arrival counts. (This was seen at EWR also, 
where the average number of cancellations per 
thunderstorm or IMC even! was more than 26 flights.) 
An airline does not cancel a flight because the demand 
it will generate might cttu»o delays. Airlines cancel 
flights because they anticipate —or are already 
experiencing -ooiily disruptions from other causes. 
Ah hough high peak, dumaod usually in creases peak 
delay, high daily demand is negatively correlated with 
high average daily delay su DFW. 

Figure 16 summands the effect of weather oo airborne 
delay ad DFW in 1997. For this figure we initially 
divided the days Into four categories; with and without 
thuadentturm* and with and without periods of IMC. 
We found that solid VMC days with thunderstorms had 
mean delays almost as small a:> VMC days without 
thunderstorms, likely because the stono.s were far from 
the airport and good visibility at the airport helped cleat 
any queues that occurred from fiow disruptions. 
Consequently it is not necessary to distinguish between 
the two types of VMC days, and figure 16 summarizes 
CODAS airborne delay tor only three weather 
combinations. 

An shown in pari a;, VMC days had the smallest 
average CODAS airborne deiay <1.9 minutes per 
oi rival). Days with IMC and no thunderstorms 

averaged 2.9 minutes of delay per arrival. Days with 
thuiKlcrstorms plus IMC averaged 6.1 minutes of delay 
by ixr arrival, more than double that of atonn-free IMC 
day*. 

Part b) shows the number of days at DFW in 1997 that 
experienced each weather category. 237 days wete 
s<>:id VMC. 79 days had one or more houri» of IMC but 
no thunderstorm activity within 5b NMI of' the airport. 
And 49 days had one or more hour* o* IMC plu$ 
thunderstorms withm 50 NMI of *he airport. 



Figure 16. CODAS airborne delay statistics for three 
weather conditions ~ DFW 1997* 

Part c) gives the resulting cumulative annual delay for 
each of the three weather conditions. Because the many 
small VMC delays occurred regularly during daily 
arrival rushes, they contributed 46% of the annual total. 
T he 79 IMC days without thunderstorms contributed 
24% of the annual total. The 49 days that had *K>th 
thunderstorms and periods of IMC contributed the 
remaining 39%, which was the second largest 
cumulative annual delay. These 49 days also included 
5 of the 6 ground hold days for flights into DFW' in 
1997. 

The predominance of c<u)g«sticm- generated VMC delay 
would imply that the FAST design does indeed focus 
on the principal queuing problem. However, analysis 
of radar data at DFW has shown that a significant 
fraction of aircraft currently land in. VMC with lea than 
minimum radar separation }* FAST cannot reduce 
separations below radar separation standards. Thus, in 
VMC periods, FAST must focus principally :>o 
planning, sequencing, and runway balancing to 
reducing large gaps in the arrival flow. 

FAST will be most effective in calm IMC period* when 
u can focus on achieving minimum separation 
standard* white also reducing arrival gaps. 

Storms can temporarily disrupt the operation ct PAST 
because it* algorithms currently depend on the use of 
standard arrival routes. FAST Is thus not able to 
achieve any of its goal* during periods when storms 
force aircraft to depart from standard routes. From 
analysis of the ITWS data log* we estimate that storms 
may have blocked arrival routes for a* much a* 50% of 
the day on the 49 thunderstorm day*: in 1997 
Fliininating these disrupted period from the modeled 
benefits for DFW would reduce the dollar saving! from 
FAST by about 7 percent 
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On the other hand, after a thunderstorm has passed, it 
often leaves a residual queue of approaching aircraft 
capable of landing at a rate close to the maximum IMC 
capacity of the airport These are the conditions for 
which small incremental increases in capacity from 
FAST can produce very large reductions in delay. 
Consider for example, a situation in which the storm 
has left scattered aircraft that now approach from many 
directions at random intervals, but with an average rate 
that is 95% of the arrival capacity of die airport Steady 
state queuing theory indicates that delay accumulates at 
a rate of 18 hours of delay per hour. But increasing the 
capacity of the airport by 5 % will reduce the rate of 
delay accumulation by more than a factor of 2. The 
effectiveness of FAST is magnified significantly. As a 
result, the small incremental delay caused by its 
absence during the storm may often be offset by the 
delay saved with its aid after the storm. 

CQMOiySllPWS 

Two prior FAST benefit models differ in their details 
and their dollar benefits for some individual airports, 
but predict comparable overall dollar savings. LMI 
modeled demand, capacity, and queuing delay more 
accurately, but generated FAST benefit estimates for 
only 10 airports. A simple steady-state model validates 
the LMI ranking results and makes it possible to extend 
the LMI dollar savings estimates to other airports. 
FAST benefits will accrue in all weather conditions. 
However, because en route and terminal airspace 
congestion causes queuing delay every day during 
arrival rushes, VMC days will contribute most of the 
annual delay benefit 

Benefit analyses based on the use of tighter inter-anival 
spacing to reduce queuing delay tend to overestimate 
FAST delay benefits in VMC when FAST can only cut 
delay cost by reducing large gaps. Benefit analyses 
also tend to overestimate FAST delay benefits during 
those periods when normal arrival routes are totally 
disrupted by storms. 

On the other side of the benefit ledger, the models 
underestimate delay benefits by ignoring reductions in 
downstream delay and ignoring periods following 
storms when large queues must be cleared. We 
conclude that these errors roughly balance each other so 
that the LMI annual savings estimates are reasonable, 
provided FAST can indeed reduce inter-arrival time by 
15 seconds at all runways. This spacing reduction, 
which is equivalent to a capacity increase of 3 to 4 
aircraft per hour, remains to be validated by analysis of 
radar data from operational tests. 
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ABSTRACT 

The National Aeronautics and Space Administration 
(NASA), working with the Federal Aviation 
Administration (FAA), is developing a suite of decision 
support tools, called the Center/TRACON Automation 
System (CTAS). CTAS tools such as the Traffic 
Management Advisor (TMA) and Final Approach 
Spacing Tool (FAST) are designed to increase the 
efficiency of the air traffic flow into and through 
Terminal airspace. A core capability of CTAS is the 
Trajectory Synthesis (TS) software for accurately 
predicting an aircraft's trajectory. In order to compute 
these trajectories, TS needs an efficient access 
mechanism for obtaining the most up-to-date and 
accurate winds. 

The current CTAS weather access mechanism suffers 
from several major drawbacks. 1 First, the mechanism 
can only handle a winds at a single resolution (presently 
40-80 km). This prevents CTAS from taking advantage 
of high resolution wind from sources such as the 
Integrated Terminal Weather System (ITWS). Second, 
the present weather access mechanism is memory 
intensive and does not extend well to higher grid 
resolutions. This potentially limits CTAS in taking 
advantage of improvements in wind resolution from 
sources such as the Rapid Update Cycle (RUC). Third, 
the present method is processing intensive and limits 
the ability of CTAS to handle higher traffic loads. This 
potentially could impact the ability of new tools such as 
Direct-To and Multi-Center TMA (McTMA) to deal 
with increased traffic loads associated with adjacent 
Centers. 

In response to these challenges, M.I.T. Lincoln 
Laboratory has developed a new CTAS weather 
distribution (WxDist) system. There are two key 
elements to the new approach. First, the single wind 
grid is replaced with a set of nested grids for the 
TRACON, Center and Adjacent Center airspaces. Each 


and the grids are updated independently of each other. 
The second key element is replacement of the present 
interpolation scheme with a nearest-neighbor value 
approach. Previous studies have shown that this 
nearest-neighbor method does not degrade trajectory 
accuracy for die grid sizes under consideration. 63 

The new software design replaces the current 
implementation, known as the Weather Data Processing 
Daemon (WDPD), with a new approach. The Weather 
Server (WxServer) sends the weather grids to a 
Weather Client (WxCiient) residing on each CTAS 
workstation running TS or PGUI (Planview Graphical 
User Interface) processes. The present point-to-point 
weather file distribution is replaced in the new scheme 
with a reliable multi-cast mechanism. This new 
distribution mechanism combined with data 
compression techniques greatly reduces network traffic 
compared to the present method. Other new processes 
combine RUC and ITWS data in a fail-soft manner to 
generate the multiple grids. The nearest-neighbor 
access method also substantially speeds up weather 
access. In combination with other improvements, the 
winds access speed is more than doubled over the 
original implementation 

DESIGN CONSIDERATIONS 

ABPOlth 

The new weather distribution design relies on replacing 
the current single wind grid with a set of nested wind 
grids and on replacing the current interpolation method 
with a nearest-neighbor retrieval method. These 
concepts will now be discussed further. 

The nested grid approach is presented conceptually in 
Figure I. The nested grids are defined for the 
TRACON, ARTCC and Multi-Center airspace. The 


* Copyright 02001 by M.I.T. Published by the American Institute of Aeronautics and Astronautics, Inc., with permission. 

’ This work was performed for the National Aeronautics arid Space Administration under Air Force Contract No. F 19628-00*0-0002. Opinions, 
interpretations, conclusions, and recommendations are those of the author and are not necessarily endorsed by NASA. 

: James Murphy is now with NASA Ames Research Center, M$-2 10-6, Moffett Field, CA. 
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nomirutl tpntiai resolution of the grids Is » nm ( - 2 km) chucked «i determine ' Met) of the nested grids should 

fur (he TRACON grid, 5 nrn {-10 km) tot the ARTCC be used. As an example, imagine an aircraft 

grid and 20 nrn f-40 km) for the Mufti-Cenler grid. It approaching tlie ARTCC from an adjacent Censer. The 

should fee timed that the grids are all aligned. That is, Mulri-Cenier grid is sized to include the furthest aircraft 

the grid {towns for ail three grids t»n be thought of as ui an adjacent Center for which a trajectory needs to be 

being placed an a uniform 1 tmi grid. generated (e.g.. to compute the Center boundary 

crossing time). The resolution of this grid it the same 
Each grid b rectangular in shape and sized in such a as the present single-resolution winds grid. As the 

way to encompass the region of interest. That is, the aircraft comes closer to the ARTCC, the aircraft's 

TRACON grul is sired to he encompass the TRACON position is checked for each retrieval to determine the 

regmn phis a buffet region around it. Likewise, the appropriate nested grid (note, sines the grids ate 

AR 1 CC grid encompasses the A.R1CC plus a buffer rectangular in shape, this check is inexpensive 

ami the Multi-Center grid extends out into adjacent computationally). When the aircraft crosses from the 

Centers a sufficient distance to allow boundary Multi-Center grid to the ARTCC grid, the wind 

crossings to be scheduled. retrievals ate then made from the ARTCC grid. 

Similarly, when the aircraft enters the TRACON grid, 
These nested grids tire used in the following way. Por die wind retrievals are then be made from that grid, 

every wind retrieval, the position of the aircraft is 


Multi-Center Grid 



Figure L. Nested wind grids for new weather distribution stheme. 
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Figure 2 shows tfw data sources for the three weather 
grub. Rapid Update Cycle \RUC) 40 kin winds and 
Integrated Terminal Weather System UTWS) 10 km 
and 2 k m winds As shown in the figure, the Multi- 
Outer grid is generated Tram RUC wind* alone. the 
ARTCC gnd is generated from RUC and ITWS 10 km 
winds, and the I'RACON grid is generated front all 
three sources. Moreover, ail thru* grids can be 
generated from RUC winds alone it ih* ITWS winds 
are not available. Finally. there can be alternate sources 
for the RUC data (nut shown;, including multiple RUC 
feeds and Tila model data. 

It should be noted that the domains of the weather 
source* and the nested grids arc independent. 
Like wine, the update rates of the weather sources ate 
independent of each other. When an update for u given 
weather source is received, the appropriate portion* of 
the affected nested grids are updated. 


In the previous section, if was shown how the nested 
grid approach allows the use of multiple wind grids 
updated from various weather i-otircet: with diffe rent 
npativi resolutions and update rates. However, a 
complication arises doc tt> the greatly increased number 
of grid pomts that need w be transmitted. In the old 
method, an interpolation scheme was used to allow the 
relatively coarse wind grids to be accessed rapidly. 
However, the memory requirements for the old method 
were vecy high on a per grid point basis. Previous 
studies dclemtined that the memory tequitetnents for 
the interpolation method made it infetuiiblc Tor 
extension to the nested grid approach. 4 ’* Accordingly a 
nmv winds retrieval method w as prop osed ns .shown in 
figure 3 


Sources 




Resolution: 20 x 20 am x 1 000’ 
Coverage; l<$00 x 1000 nm *40,000’ 


ART'CC: 

Resolution: 5 x 5 nro x i000’ 
Coverage: 800 x 500 xm x 40,000’ 




Resolution; i x 1 nm x 500* 
Coverage: 100 x 100 nm x 1 8,000* 


Figure i. Data source, for weather grids. 
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figure 3. Illustration of replacing Interpolation with neighbor retrieval method* 


In the new method, the interpolation step is bypassed in 
favor of a nearest-nelghbor retrieval ache me. The 
primary motivation for the nearest-neighbor approach is 
to make the memory requirements of the nested grids 
feasible, but. it also has the advantage of speeding op 
accost time and reducing processing roguirementt, A$ 
reported previously, the weather data taxes* speed more 
than doubled with the new method {note: include; the 
e fleet of eliminating geometric altitude, which is being 
incorporated into the present system;/ Figure 4 
illustrates these weather data access speed 
bnp.TJYcmenls. 

in order to employ the nearest-neighbor technique, it 
was necessary to verify amt trajectory accuracy would 
noi be impacted. As repotted in, teats ware run 
comparing tnycctcrics computed uoing the interpolation 
and nearest-neighbor methods/' The compariiion was 
tun tor two coses' a; Meter tlx to threshold using ITWB 
2 km winds and b) Coordination fix to meter fix using 
KUC interpolated ?.o iO km A 1TWS 10 km wind). 

It waii found that the use of the nearest -ntrigttboi 
tut? bod produced a one second RMS difference in 
trajectories for the first case and a four second RMS 
difference tor th« second case. These differences ore 
negligible for tho trajectories examined. An example of 
these results tor ITWS 3 km dma is ahown in Figure V 
Note the effect of using the nearest-neighbor method 
tot tlie Mult? ‘Center grid was noL examined. however. 


the effect is assumed to be fosignUkiKt given the iarge 
distances involved and low update rate <'f the adjacent 
Center traffic data (e.g., 3 minutes via ETMS), A 
recent study showed iitc cftccs of using nearest- 
neighbor access for 40 am winds for trajectories in the 
Center was » ten sscoad RMS difference vs. 
interpolation. 1 

SOFTWARE DESIGN 

ttmrto 

A block diagram of the Weather Duiixibutioa (WxDisi) 
software is shown m f igure d. The key modules are the 
Weather Server (WxBervet) and Y.'eather Client 
(WxClieni) modules. The WxSctver module provides 
the weather data to multiple WxChent modules. There 
in on WxServer module for u given CTAS installation, 
and there L> one VxClitnt module for each workstation 
employing one or more IS or PGUi pfcK<??.ses. The 
external weather sources are converted into GRIB 
(Gridded Binary) tiles and divided up into minor grid:? 
as described in the previous section. ITui miuoi grids 
arc transmitted to the WxOieot processes via the 
reliable multi'Uwt protocol ovtr the local area network 
(LAN). The WxOiem processes provide the weather 
data to the application processes via a shared memory 
interface. Die Weather Library fWx Library) s< ccsscu 
the shared memory and provides die interface to the 
weather users. 
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Figure Is Weather data access speed improvement 



Figure 5. Comparison oi l r« s fin mtrrpotMed v* ntjnxt neighbor method*. 
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n«»n> d. Weather distribution system block diagram* 


Figure 7 shows a more detailed view of t:hc CTAS 
weather distribution modules A** shown in the figure, 
the I HVS data is processed by die iTWS Connection 
Module into GRfB files and passed to the WxScrvc*. 
The WxScrvet module ah© receives RUC tiles in GRIB 
formal from the existing. WDAD (Weather Data 
Acquisition Daemon) process. Note that ffca processes 
for acquiring the iTWB and RUG data reside outside 
the firewall to isolate the weather souues front the 
CTAS aysiem. Thu CTAS Weather Communicatioxi 
module implements me reliable multi-cast protocol tor 
tmosfetrmg the minor grids over the LAN. Finally, the 
WxLibrary module provides the interface between tlw 
weather dtaibution system and rite *q>pi unions 

Tahiti I lists the major modules in the owv weather 
doanbutioo uystem. The functionality ol each module 
is summarized in the tabte. These modules will now be 
Me fly described. 

EMrj?Matt«lato.,8;to» a te 

The ITWS Common;; ation Module fitw.i rtdc) obtain* 
the Terminal Winds data feed from the ITWS testbed 
via stream connections and converts the data from 
QirUwutn (CAR) format to GRID (Gudderi Dinury) : 
format. It establishes .socket connections to tm/cmit the 
ennvuned JTWS Terminal Winds data to multiple local 
and/or remote Weather Server*. 


Thrt Weather Server (WxServer) module sets op the 
aiie-adaptabfe wind grids from u WxSeivcr 
configuration file*. 1 he Hie d*ilne» tho resolution, 
spatial extent and data sources for each grid. It also 
defines the backup strategy for generating ih»$< gridj 
from alternate sources. 

T he Wx Server accepts Terminal Winds data from the 
JTW5 CommumuAlon^ module via socket connections 
and RDC data from this existing Weather Data 
Acquisition Daemon (WDAD) process via file transfer. 
The WxServer automatically selects between the 
alternate weather sources to determine the best 
available data. It rteftt translates ihu mpue data to the 
NAS coordinate system sud comprise ses the dal* for 
U^namiVton usinz the GRIB compression algorithm. It 
then transmits the data to each Weather Client m.lng the 
reliable muhica.r protocol. 

The demils of the automatic source selection lo^fc is 
shown in Figure #, In this example, the ITWS 2 km 
winds data feud is interrupt'd at time T*. A tier a 
iimeout period (nominally 6 minutes), this weather 
tountt *s dec lared unavailable. Howessrr, the la?t data 
received com in net to ht? bsod until the nominal ITWS 
10 km update urns. 



If the ITWS 10 km update is received, the ITWS 2 km 
grids are now ge*>erated from the 10 km winds data. If 
the 10 km data is not available, then this source Is 
declared unavailable and the most recent data continues 
to be used until the nominal RUC update time. When 


the RUC update is received, then the ITWS 2 km and 
10 km grids are now generated from the RUC data. If 
the ITWS data later becomes available, then the system 
returns to using that data to generate the 2 km and 10 
km wind grids. 



Figure 7. CTAS weather distribution modules. 


Table 1. Major Software Modules 


NAME 

FUNCTIONALITY 1 

ITWS Connection Module 

• Connection to ITWS for winds data 

• Data conversion to GRIB format 

• Connection to multiple WxServers 

WxServer Module 

• Site adaptable weather grids 

• File connection for ITWS data 

• Conversion to NAS coordinates 

• Automatic weather source selection 

• Weather data compression 

Communications Protocol 

• Automatic conncction/reconnection 

* Reliable multicast to transmit compressed weather data 

WxClient Module 

• Read multicast weather data 

• Update shared memory buffers 

• Switch buffers on command 

WxLibrary Module 

9 When user asks for weather products: 

1 . Determine appropriate grid 

2. Select nearest-neighbor value 

3 . Read weather product 

• Switch memory page on command 

• Identify current weather sources 
Note: grid structures transparent to users. 
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Figure 8. Tixnrftoess check logic iltmd ration. 


Figure 9 summarizes the Wx&erver/WxClient 
communications protocol. The VVxServer process 
dehvers the nested ^rid date to cue or more WxClient 
pt<*oeii?.c$ usmg. a reliable multicast protocol. This 
protocol makes us* of a TCP/IP based socket 
connection between the Wx Stiver and each Wx&ient 
as a control channel. i\ uses IP Multicast as the dam 
channel from the Wx Server to all \V* Clients. 

The use of TCP/IP tot the control channel provides 
reliability and ensures that no control message will Ikt 
lost or delivered out oi* serpens*. this guarantees the 
liability of the shared data view within the 
W;<Set ver/WxCT tents group, The u^e of IP Multicast 
for the data channel allows us to transmit the data once 
for ail WxCiisnta legatdless of how many axe in use on 
the CI A S system. Tins minimises the network load 
associated with weather data and allows the 
transmission of higher spatial and temporal resolution 
data seta. 

Mite 

The Weather Client fWxCiicnt) rcada the weather grids 
via the reliable omit [cast protocol and updates the 


appropriate atea of Stated memory . The details of the 
shared memory interface with the weather users is 
discussed below. 


Shared memory is used by the W.\Clteni to make the 
weather data available the the Wx User, Two buffers 
are used for each data grid, one contains the current 
data and is; available tor reading, the other is being 
written to with iww data. The form of each buffer is 
fixed wish a standard header followed by a three 
dimensional array of product structures. 

Ihe doable buffered approach used for shared memory 
is illustrated jrt Figure >0. As shown in the figure, new 
data is written into the wiitc page while the users stucco 
the mad page. When the updae Is complete, the pages* 
are swapped. There is a signal handling scheme 
employed to ensure that all the users have the current 
information before the swap Is cai/ied out. Another 
implementation detail i:> then the contents of the new- 
read page nerd to be copied back to the oid read page 
prior to allowing further updates to occur to the new 
write page. 
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Snpport for Multiple Weather Users 

Each WxCUent is capable of supporting multiple 
weather users. Recall that there is one WxClient for 
each workstation that supports all the weather users 
residing on that physical piece of hardware. As 
described in the previous section, a double buffer 
scheme is used to allow the weather users to read the 
current weather data from one buffer while the other 
buffer is being filled with the next weather data to be 
used. An important consideration is to ensure that all 
weather users switch from one buffer to die other buffer 
in concert This synchronization is carried out using a 
semaphore mechanism. 

The Weather Library API provides the necessary 
interface to the weather server for all weather using 
applications. The intent is to allow the weather using 
application to request one or more weather products at 
one or more locations in a single request Units used 
for both the locations and the weather products are felt 
to be those most suitable to the using application. 

ALGORITHMIC ISSUES 

There are three issues which have been identified as 
potentially requiring changes to T$ in order to 
accommodate the new weather distribution system. 
These issues include: 


be very large. The TS software copes with large 
gradients by limiting the vertical wind gradient to 10 
knots. The zero gradient simply means that the effect of 
wind is ignored. In neither case does TS fail in its 
trajectory calculations. 

One specific instance where the vertical gradient Is used 
is in die en route portion of the trajectory: specifically 
the Constant CAS and Constant Mach segments. As 
pan of the TS algorithm, a system of differential 
equations is solved using a discrete step method called 
Rimge-Kutta. This applies in particular calculations of 
True Airspeed, vertical speed and ground speed. 

The principal equation in TS in solving for V, is (Eq 1): 



T-D 

m 


m Cu 


( 1 ) 


where y a is the aerodynamic flight path angle, T is 
thrust, D is drag force, m is aircraft mass and w its 
weight, V w is the wind speed and 6 r w is the relative 
wind angle. 

For constant Mach or constant CAS, the above 
differential equation reduces to an algebraic expression. 
This is apparent if one considers the fundamental 
relationship (Eq. 2): 


• gradient computation, 

• temperature interpolation and 

• capture condition completion. 

These issues will now be discussed. 

The first TS weather use issue is the calculation of the 
vertical wind gradient This gradient was originally 
calculated by accessing wind data at two slightly 
different altitudes h and h+Ah, taking the difference and 
dividing by Ah to obtain a discrete approximation of the 
vertical wind gradient 

In the original system, a value of 50* was used for Ah, 
which works properly when interpolation is used. 
However, when nearest-neighbor retrieval is used with 
this small value, the coordinates h and h+Ah usually 
round off to the same grid point and produce the same 
wind values. The resulting gradient is then zero. 


V, ma(h)M(y CAS ,h) 

where a(h) is the speed of sound (as a function of 
altitude) and M(„.) is the Mach number (as a function 
of Fa. and altitude). It is clear that when Mach or CAS 
is constant, that the true airspeed is a simple function of 
altitude. 

One can then write (Eq. 3): 

dV, _ dh dV, (3) 

dt " dt‘ dh 

where (Eq. 4) 



Making the approximation that the wind velocity over 
short x or y distances is constant (i.e. only dependent on 
h) then (Eq. S): 


In an exceptional case, the second value would round to 
die next altitude level and the resulting gradient would 
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jIT.casOJ - (5) 

This then leads to the expression for the aerodynamic 
glide slope (Eq. 6): 


Ya dV ll (6 > 

ah ah 

Note that the wind gradient shows up in the 
denominator. Thus, die reason for limiting the gradient 
value to a maximum value is to avoid a singularity 
condition for the aerodynamic glide slope and 
eventually a sign reversal. This would certainly make 
TS fail, since in a descent phase it then could not satisfy 
the required boundary conditions. It can also be seen 
that if the wind gradient comes out to be zero, the only 
effect is that the glide slope used in the trajectory 
prediction is a slightly small. 

But the same formula uses weather in other ways: the 
altitude derivative of the True Airspeed uses 
temperature readings at two different altitudes. This 
would lead to problems also if the chosen altitude 
difference is a small value. 

The solution for both die wind gradient and the True 
Airspeed gradient is to force the altitude increment to 
be equal to the weather grid vertical spacing. In the 
case of nested weather grids that increment could be 
different for the different nested grids. This introduces 
some complications in the event that the gradient 
computation crosses the nested grid boundaries. For 
this reason, it would be preferable to introduce an 
explicit call for die gradient at a particular point. This 
would allow these complications to be isolated from the 
weather user. 

Temperature Interpolation 

The second TS weather use issue involves temperature 
interpolation. One important TS function is to meet 
capture conditions, such as matching cruise and descent 
segments to identify the top of descent (TOD). In 
certain cases involving iteration to meet capture 
conditions, the TS was found to fail when nearest- 
neighbor retrieval was used. This is because TS uses 
small changes in temperature with altitude to drive the 
solution in the correct direction. When the nearest- 
neighbor value is used, the TS could possibly to 
converge because the same temperature value is always 
retrieved. 


This problem was observed for the ARTCC grid when a 
grid spacing in altitude of 2000* was used. However, it 
was found when the altitude grid spacing was reduced 
to 1000*, no TS failures were observed. As previously 
noted, the grid spacing in altitude can be traded against 
the horizontal grid spacing without increasing the 
memory required if even greater vertical resolution is 
needed (e.g., the vertical grid spacing could be 
decreased from 1000' to 2S0* if the horizontal spacing 
was increased from 10 am to 20 nm). 

Another considered was a fix to always linearly 
interpolate the temperature values between the nearest- 
neighbor values above and below the current altitude. 
This approach incurs a minor performance penalty due 
to the need to retrieve two temperature values (instead 
of one) and to perform a simple interpolation. 
However, this approach has the virtue of keeping the 
API unchanged. It also improves the accuracy of the 
temperature values, which feature a strong dependence 
on altitude. Howeva , in practice it has been found that 
the decrease in vertical grid spacing proved sufficient 
and this fix was not implemented. 

Capture. GmMm 

A problem was found in the way that TS completes the 
capture condition iteration. . When TS iterates to point 
where the altitude is within the capture limits, h stops 
and returns the capture altitude. For example, the 
desired capture altitude might be 25,000' and the actual 
capture altitude might be 25,010'. Originally, TS did 
not recompute the derived variables (CAS, ground 
speed, etc) for the desired capture altitude but returned 
the derived variables for the actual capture altitude 
instead. 

This behavior caused problems In the TS computations 
due to the nearest-neighbor retrieval for temperature. 
Fore example, if the gridded winds layers had 2000* 
vertical spacing, then the temperature might be 433 °R 
at 24,000* and 426 °K at 26,000*. The nearest-neighbor 
temperature value would therefore be 433 *R at 25,000* 
and 426 °R at 25,010*. This would create a 
discontinuity in the temperature values between flight 
segments that might cause TS to fail. 

A fix was implemented is to force TS to recompute the 
derived values at the end of the iteration for the desired 
(instead of actual) capture altitude. This fix has been 
accepted for incorporation into the CTAS baseline 
software. It should be noted that the proposed 
temperature interpolation fix would also address this 
particular problem. 
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TESTING PROCEDURES 

Tests were carried out on the prototype implementation 
to validated functionality and measure performance. 
Note: the results presented here should be considered 
preliminary in nature and subject to further refinement. 

Foncttonalltv Testing 

The functionality tests were carried out in several steps. 
The first step was to perform a regression test using 
RUC data only. The second step was to add ITWS 
winds and validate correct insertion into the wind grids. 
The third step was to validate that the wind grids 
continue to be properly generated when the ITWS 
winds art transiently added and removed. The fourth 
step was to verify proper TS operation in the presence 
of wind field discontinuities. The fifth step was to 
quantify the difference in ETA (Estimated Time of 
Arrival) values with die addition of ITWS winds. 

For testing purposes, a version of the new system was 
created which returns the linearly interpolated weather 
value instead of the nearest-neighbor value. This 
version (WxDist Interpolated) is not intended for 
operational use (since it runs more slowly than the 
nearest-neighbor version) but allows direct comparison 
between the weather values returned by the new vs. old 
systems from the getWeatherValue function. 

The CTAS software was also instrumented to generate 
ETA logs, ETA log summaries, track logs and 
getWeatherValue logs. The standard output is to 
produce the ETA log, ETA log summary and track log 
(if radar track data is used). If verbose output is 
selected, then the getWeatherValue log is also produced 
(generally limited to short runs due to the large volume 
of output generated). 

A capability to generate synthetic RUC data sets was 
also implemented to assist in regression testing. A 
utility program was written which allows synthetic 
RUC data to be generated for various test conditions. 
For example, one file was created with uniform wind 
values and temperature values that increased, linearly 
with pressure level, and a second file was created with 
uniform temperature at all levels and U &V values that 
increased linearly with RUC X & Y coordinates, 
respectively. These files were used to validate the 
RUC-to-CTAS coordinate transformation in the vertical 
and horizontal dimensions, respectively. 

Additional utility programs & Unix scripts were written 
for examining the input RUC data and processing the 
output test data. Unix scripts were also written to 
simplify making test runs and to automatically save the 


output logs. Processing and examination of the test 
data was primarily done using IDL and Excel. 

Regression Testing 

Regression testing was carried out to ensure that the 
new weather distribution system preserves the CTAS 
functionality. In particular, it is necessary to ensure 
that trajectories are correctly generated with the new vs. 
old systems. In order to do this, it is necessary to use 
RUC data only, since the old system cannot ingest 
ITWS winds data. 

RUC winds Ingestion 

The first regression test was to verify that the RUC data 
is properly ingested into the new system. In order to do 
this, actual and synthetic RUC riles were input to the 
WDPD, WxDist Gridded (nearest-neighbor) and 
WxDist Interpolated versions of CTAS. The outputs of 
the three versions were then compared and any 
inconsistencies diagnosed. In the process of carrying 
out the regression testing, several problems were found 
in the WDPD processing which were diagnosed and 
fixed. These changes are being evaluated for 
incorporation into the current CTAS baseline software 
but will not be further discussed here. 

Examples of this comparison are shown in Figures 11 
and 12. For these examples, the getWeatherValue calls 
were logged using the verbose option for a single 
aircraft trajectory. Figure 11 shows the retrieved 
temperature vs. altitude and Figure 12 shows the 
retrieved U wind vs. altitude. As seen in the figures, the 
WxDist gridded (nearest-neighbor) retrievals exhibit 
the expected staircase behavior whereas the WDPD and 
WxDist interpolated values vary smoothly. (Note: there 
is a small anomaly in the WDPD results below 5000' 
which is currently being investigated.) 

FMwnpairfoon 

The second regression test was to compare trajectories 
generated using the WDPD and WxDist Gridded 
weather distribution schemes. For this test, 489 flight 
plan trajectories were compared for 158 aircraft using 
DFWF RUC weather data from November 22, 2000 at 
1400Z. Figure 13 plots the ETA difference between the 
new and old systems for meter fix to threshold 
trajectories. The worst case differences range from -7 $ 
to + 15 s. The mean ETA difference was 0.9 s (0.12%) 
and the RMS ETA difference was 3.3 s (0.44%). Also 
shown is the comparison for the WxDist Interpolated 
version which produced identical results to the WxDist 
Gridded version. 
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Figure 13, flight Flan ItA difference between WDFf> and WxDivr for m titer 0.x to threshold trajectories. 


Figure 14 shows the ETA difference tot the same tci of 
aircraft for trajectories from the coordination fix to the 
toeicr fix. The worst* case FT A differences range from 
•45 s to +31 The mean KTA difference was *1.0 s 
(0.06%) nod the RMS difference was 5.6 :» (0.3>%). 
Generally, the values matched very closely evtih the. 
exception of a few luige diffcicticcrc which are currently 
under investigation. Also shown is tire comparison for 
the WxDtst Interpolated version which were essentially 
the same as the WxOitt Gridded results except that 
there was one fewe* large excursion. 

These error differences arc slightly larger than 
predic ted by reference (6j. That study predicted ETA 
ddTcteiKea of I. s RMS for ft* meter fix to threshold 
erase and 4 RMS fot die coordination fix to meter fix 
cate based on the use of ngtreat-netgbbor vs. 
interpolated weather value retrieval, However, the 
ETA differences between the WxDlst Ciridded and 
WxDitt Interpolated versions arc in fact much smaller 
that predicted by the study, The observed variation b 
therefore likely to be dun to a;; yet undiagnosed 
differences between WDH) and Wxiaisi However, die 
test results show that these differ cnees have been 
reduced to very small values, 

A comfwuison between the new and old systems was 
also made fur aircraft trajectories troto re-lar tracks. 
Figure 15 shows tlie Time To Fly (TIT) horn 426 rodnr 


track trajectories for a » aircraft la ruling on runway t?|, 
for this case, it should he noted that the WxDist ETA 
values were j&AJ&gj for the gridded and interpolated 
versions. The worst-case differences in TIT (ETA- 
current time) were *7 s to t5 seconds. The mean 
difference was I.) a (0.7.3%) and the RM.S difference 
was 1.5 s (QJ$%) Additional regression testing is it) 
progress, hut these result* suggest that the new and old 
systems arc working nearly identically fot radar tracks. 

11.V8? Wtedii VitMittofi 

Jlic next reeling step was to validate that ITWS HI km 
and 2 km winds are inserted correctly into the: wind 
grids. For this tost, the wind grids were defined as 
shown io Figure 16. For these test*, separate major 
grids were used for the TEWS 10 km »ml 2 km data. 
These major grids were set to correspond to the 
maximum extent of the XTWS wind fields. 

in order to readily verity ptoper insertion of the fTWS 
winds, a synthetic GRIB data tool was ueed to generate 
dummy HTVS wind files. Two fdcs were created, one 
for ITWS 10 kin with uniform winds due (Taut sad the 
other for 1TWS 2 km wbh uniform winds due North. 
These ITVv'S wind files were then processed by the 
Weather Server with tie* RUC wind file for 25 August 
2«KKi af t#00/ for 5300' (350 mb) hi produce the 
multiple wind £iids. 
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Figure 16. Major grid definitions for XTWS winds integration test 


The wind field values were retrieved on a grid 
approximating the 40 km RUC grid. These results are 
as shown in Figure 17. As seen in the figure, the 
insertion of die dummy ITWS wind data is clearly 
demonstrated. 

Transient ITWS Winds Availability Testing 

The third functionality test was to verify proper 
operation with transit". ITWS winds availability, litis 
test was carried out by interrupting the ITWS winds 
availability during a nonnal run and verifying that all 
the major grids continued to be generated from the 
RUC winds only. 

Wind Field Discontinuity Testing 

The fourth functionality test was to verify proper TS 
operation in the face of wind field discontinuities 
between the RUC and ITWS data. In order to rule out 
this possibility, CT AS was run with the dummy ITWS 
files shown in Figure 17. The results were examined 
and showed no evidence of TS failures in the face of 
woist-case discontinuities. 

Effect on ETAs of Including ITWS Winds 

The fifth functionality test was to quantify the effect on 
ETAs of including ITWS winds, as shown in Figure IS. 
To carry out this test, CTAS was run with RUC-only 


winds vs. RUC plus ITWS winds. Earlier studies 
showed considerable variation in meter fix to threshold 
ETAs as a function of time should be observed with the 
inclusion of ITWS 2 km winds updated every five 
minutes. 5,7 

Figure 18 summarizes the result of computing 489 
pFAST (Passive Final Approach Spacing Tool) 
trajectories from 158 arrival flight plans under two 
conditions: 1) RUC winds only with interpolated winds 
(old system) and 2) RUC + ITWS winds with gridded 
winds (new system). These results were computed for 
DFW on 11/22700 from I500Z to 1700Z with RUC 
winds updated hourly, ITWS 10 km winds updated 
every 30 minutes and ITWS 2 km updated every 5 
minutes. 

As seen in the figure, the ETA values for the RUC only 
winds increase by 3 seconds over the two hour period. 
By contrast, the ETAs for the RUC + ITWS winds 
change by 16 seconds over this time period. These 
results are consistent with the earlier studies. 

EttfopnfiiittXgtiM 

Performance testing was carried out to assess the 
processing speed and memory requirements for the new 
vs. old systems. Table 2 shows examples of 
preliminary results for various grid sues. 


46 




Vi3 


AIAA CuUlant-t, Navigation A Control Conference, Montreal, Quebec, W August 2WI 


( 




Figure IH. pfACT ETA comparison utntmar) for MJC/lnterpolaM vs. RUC-HmSAirnltict) winds. 


4 


A1 AA Guidance, Navigation & Control Conference, Montreal, Quebec, 6*9 August 2001 


Table 2 . Example performance measurement results (preliminary). 


Process 

Total Data 

Points, 

XxYxZ 

CPU Use. 
seconds 

Memory 

Network Load, Mbyte 

Use, 

Mbyte 

At the 
Hour 

At the 

Half 

Hour 

At 5 
Minutes 

WDPD* 

30x16x55 

11.73 

24 

30** 




15x8x55 

4.86 

18 

7** 




9x5x55 

3.37 

18 

2** 



WxServer* 

30x30x50 

10.98*** 

84*** 

0.7*** 

0.5*** 

0.2*“* 


24 x 24 x 50 ] 







70 x 70 x 70 






WxClient**** 


8.85*** 

34*** 





•270 MHz Ultra 5 - 128 MB RAM 

** Assumes 5 TS, 1 PFS, 2 PGUIs = Eight Using Applications 
*** Total for all grids 

•**•60 MHz SparcStation 20 - 64 MB RAM 


SUMMARY 

This paper described die design and implementation of 
the new CTAS weather distribution system. The 
approach relies on two key concepts. The first concept 
is the use of multiple nested wind grids for the 
TRACON, ARTCC and Adjacent Center airspaces to 
replace the single low-resolution weather grid currently 
used. This new method allows the use of higher spadal 
and temporal resolution products such as IWS Terminal 
Winds and improved RUC winds. The spatial 
resolution and update rate for each grid is tailored to the 
weather sources and user requirements. 

The second concept is the use of nearest-neighbor data 
retrieval to replace the interpolation method currently 
used. Previous analysis showed that the memory 
requirements of the present method prevent its 
extension to higher resolution weather grids. The 
nearest-neighbor method was introduced in order to 
reduce the memory requirements to feasible levels. As 
an added bonus, the nearest-neighbor method also 
yields a substantial improvement in weather data access 
speed. Previous work showing that use of the nearest- 
neighbor method should not substantially degrade 
trajectory accuracy for the grid resolutions under 
consideration was confirmed in the present study. 


The software design employs the concept of a single 
Weather Server process that provides weather data to 
multiple Weather Client processes. There is one 
Weather Server for a given CTAS site installation and 
one Weather Client for each workstation running one or 
more T$ or PGUI processes. A reliable multi-cast 
protocol is used for transmitting the weather data from 
the Weather Server to the Weather Clients. The 
weather grids are divided up into subunits (called minor 
grids) and compressed for transmission using the GRIB 
format data compression technique. Each minor grid 
has a primary weather data source and optional 
secondary weather data sources. In the event that the 
primary weather data source is not available, the minor 
grid can continue to be generated using the secondary 
weather data sources. 

The software is divided into ITWS Connection, 
Weather Server, Communications Protocol, Weather 
Client, Weather Library and Weather User modules. 
The ITWS Connection module inputs ITWS Terminal 
Winds and converts it to GRIB format. The Weather 
Server merges the RUC data (from the existing WDAD 
process) and the ITWS data to generate the nested grids 
information. The Communications Protocol module 
performs the reliable multi-cast of the nested grid data 
to the Weather Client processes. The Weather Client 
module receives the nested grid data and makes it 
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interface (API) for the user processes which selects the Lincoln Laboratory, 18 July 2000. 

appropriate nested grid for data retrieval in a 

transparent manner. The Weather User module 

represents changes to the application processes where 

needed to accommodate the new method. 

The testing procedures were also described. Results 
were provided including regression testing, 
performance measurement and software metrics. 

Finally, future work was described. 
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ABSTRACT 

Knowledge of present and future winds and 
temperature is important for air traffic operations in 
general, but is crucial for Decision Support Tools 
(DSTs) that rely heavily on accurately predicting 
trajectories of aircraft One such tool is the Center- 
TRACON Automation System (CTAS) developed by 
NASA Ames Research Center. 

The Rapid Update Cycle (RUC) system is presently the 
principal source of weather information for CTAS. U 
RUC provides weather updates on an hourly basis on a 
nationwide grid with horizontal resolution of 40 km and 
vertical resolution of 25 mb in pressure. 1 However, a 
recent study of RUC data availability showed that the 
NWS and NOAA servers are subject to frequent service 
interruptions. Over a 210 day period (4/19/00- 
11/11/00), the availability of two NOAA and one NWS 
RUC server was monitored automatically. It was found 
that 60 days (29%) had periods of one hour or more 
where at least one server was out, with the longest 
outage lasting 13 hours on 9/21/00. In addition, there 
were 9 days (4%) for which all three servers were 
simultaneously unavailable, with the longest outage 
lasting 6 hours on 5/7/00. Moreover, even longer 
outages have been experienced with the RUC servers 
over the past several years. 

RUC forecasts are provided for up to 12 hours, but 
these are not currently used in CTAS as back up 
sources (except that the 1 or 2 hour forecasts are used 
for the current winds to compensate for transmission 
delays in obtaining the RUC data). Since RUC outages 
have been experienced for longer than 12 hours, it is 
therefore necessary to back RUC up with another 
weather source providing long-range forecasts. 

This paper examines the use of the Eta model forecasts 
as a back-up weather source for CTAS. A specific 


"This work was perf o rmed for t be Federal Aviation Administration 
under Air Force Contract No. FI 9628*00*00002. 
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output of the Eta 32 km model, namely Grid 104, was 
selected for evaluation because its horizontal and 
vertical resolution, spatial extent and output parameters 
match most closely those of RUC. 4 While RUC 
forecasts for a maximum of 12 hours Into the future, 
Eta does so for up to 60 hours. In the event that a RUC 
outage would occur. Eta data could be substituted. If 
Eta data also became unavailable, the last issued 
forecasts could allow CTAS to continue to function 
properly for up to 60 hours. 

The approach used for evaluating the suitability of the 
Eta model and RUC forecasts was to compare them 
with the RUC analysis output or 0 hour forecast file, at 
the forecast time. Not surprisingly. It was found that 
the RUC model forecasts had lower wind magnitude 
errors out to 12 hours (the limit of the RUC forecasts) 
than the Eta model had. However, the wind magnitude 
error for the Eta model grew only from 9 ft/s at 12 
hours (comparable with RUC) to 11 ft/s at 48 hours. 
We therefore conclude that RUC forecasts should be 
used for outages up to 12 hours and Eta model forecasts 
should be used for outages up to 60 hours. 

MsmomoGY 

The comparison of RUC and Eta data was done for a 
typical ARTCC (Air Route Traffic Control Center), in 
this case the DFW Center (2FW) airspace, in the 
vertical dimension, three altitude layers were examined: 
7,500*, 18,400* and 30,000*. The time period over 
which the comparison was made was a period of ten 
days, starting 9 June 2000. 

Since the RUC and Eta model data is provided on 
different projection systems (Lambert Conformal vs. 
Polar Stereographic) and on different grids, the Eta data 
was first transformed onto the RUC grid using the 
following procedure. 9 Each RUC grid point was 
transformed into the Eta grid system and the 
surrounding Eta grid points determined. A linear 
Interpolation was then done using the Eta weather 
products at the eight comers of the cube surrounding 
the RUC grid point The resulting interpolated Eta 
model value was then used for the corresponding RUC 
grid point 
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Once i lie data were on ihc same grid, metric* were 
adopted for comparing weatbef data. For tcmpct*iure. 
tile metric was simply the scalar different* between fhe 
RUC analysis teroperatorc and me RUC or £ta modd 
forecast temperature. For winds, there art* two possible 
nmnes; vector difference and ni&gmtuite/dtreaion 
difference as iflusirated in Figure I . 


Difference Vector 



AV Magnitude Difference 
AO Direction Difference 


Figure L Metrics for comparing wind vectors, 

Tlxc vector difference {shown by the dashed line in 
Figure ); is often used us u xWnd error metric. 
However. thW metric provides limited insight since the 
ouemeUcn of the error vector varies greatly and iho 
magnitude of the error vector is strongly influenced by 
the difference in direction between the two original 
vectors to be compared. For this study we chose the 
mag ni tude/di rection metric tAV,A©; as being more 
physically meaningful 

Finally, a selection must he made of the weather files to 
he compared. RUC outputs a set of forecasts on an 
hourly cycle, while Ufa outputs a set of forecast!; ott a 
s'u hour cycle and while the. louecat RUC forecast 1$ 12 
hours, Bta forecasts as far as <S0 hours into die future. 
Figure 2 shows the RUC arid Eui cycles and the set of 
files generated at each update. The find file in euCh set 
»s called ihc 9 aualystfi file’ 1 f abbreviated an I or tO if for 
RUC and cO if for Etay 'Thu; file represents the beg? 
Knowledge of the weather at iho analysis time after all 
the new measurements gathered during the previous 
cyde have been incorporated. The files following thut 
are forecast: a vlxhour RUC forecast should be labeled 
fQ6 or simply .*6 for RUC and e6 for Fra. A comparison 
between the RUC analysis flit: and an Eta la bour 
forecast w-MiId then be tteorecri by rfife IS cm 


RUC Output Schedule 

1 !*«*>< S <^RMWV 



Eta Output Schedule 



Figure 2, RUC and Eta made! output schedules. 

The comparisons that were made are first between the 
RUC and Eta analysis files, denoted .rO/eO, for sdl sets 
output at the times corresponding to die Eta cycles (i,e. 
every 6 hours) for u period of 10 days. The 40 
comparisons thus made were statistically analysed 
separably and then combined for a global result in the 
end. Similarly, comparisons between RUC forecasts 
and a reference weather file, and Eta forecasts arid the 
same reference file were made for me same 40 cycle 
times aepn.tai.cly first and then the results were 
combined for a global result. The reference tries were 
chosen to be the RUC analysis tiles that were generated 
at the time to which the forecasts were ptojened. i.e. a 
12 -hour forecast was compared with the ro file 
generated 12 hours later. Although this seems to 
disadvantage Eta in a dura quality comparison, it 
scorned justified by the fact that Eta is only u backup 
system to RUC and that RUC is the system normally in 
use. Finding an independent went her source (such as 
MDCRS readings etc) other than RUC or Urn io play 
ihe rule of ground truth mold he done but wan beyond 
ihe scope of ibis study . 


52 



AIAA Outdance, Navigation <& Control Conference. Montreal, Quebec. Augu$t H 2001 


KU 


Figure 3 presenu the results for the rO/cO, to-called 
'direct* data companion The *-axi.i is the index to the 
40 Eta cycle update* and therefore the number of 
independent comparisons. The window* repaint the 
results for compnxisou of temperature, wind strength 
and wiud direction, averaged over the ZF W area and 
over the thccc selected altitudes. 
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Figure 3, Direct RUC vs. Eta model t ompar istm of 
tempe rating wind strength mid wind dim lion for 
40 hourly np<J;U< v (analysis rytlrs), 


fn order to better visualize what these results 
summarize we show in Figure 4 an overlay of the RUC 
and Eta wind fields at the altitude of 30,000ft 
OQOmhar). One observes immediately that both 
systems present the same weather pattern, and that the 
wind strengths match well, but that there is a very 
discernible difference in wind direction. This seems to 
hold true for all comparisons made, whether among 
analysis files or forecasts. Wind strength varied from 5 
It/s in some parts to 45 ft/s in other parts, yet the 
average differences, returning to Figure 3, hover around 
zero, with a StD (standard deviation; ci around <5ft/s for 
wind strength difference and a StD of about 1,25 
degrees F in temperature difference. This seems to 
indicate that the differences between the RUC or Eiu 
analysis files are small and that we are justified in 
taking either one (we choose RUC) as our reference 
when evaluating the RUC and E la forecast*. 
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Figure 4, Overlay of RUC and Eta wind vectors. 


1 he neat comparisons art between RUC forecasts and 
the KDC analysis data, generated at the time to which ia 
forecasted; and similarly between Eta forecast* and 
RUC analysis weather that riraaetiati/od later at the 
appropriate lirtMX 


Figure $ shows the results for the labour forecasts of 
the temperature: one curve for r0/rl2 and another for 
r0/c 12, with as x-axis the 40 Elm update time*. In tiic 
first window we show the mean difference and in the 
second the 5tf> of the temperature difference. In 
figure* b und 7 we show similar results for the wind 
strength snd wind direction differences. A close look 
reveals that peaks and valleys in both curves match 
quite closely. This means thru there were weather 
changes not predicted by either system, and that they 
erred in the similar ways. 
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Figure 8 shows means and StD for Use comparison 
rO/t4H t for the 4 Sftour Bia I’orec**! (no curve for RUC 
can he shown since RUC only fotccasu for up to 12 
hours). 


Figures 9,10 and 1! summarize all tltese results 
averaged over till 40 u/xJate*. These new figures contain 
the mean and StD of diHeftoces rt)/ri for 1,2,3 ,6 and 12 
hours and rO/ej for 0,6,12.18*..48 hours fot the 
parameters temperature, wind strength and wind 
direction. 
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Figure 5: Comparison of 12-hour temperature predictions fur RUC and Eta forecasts 
vs, RUC ana lysis data for 40 analysis cycles* 
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Figure 6: Comparison of 12- hour w ind strength predictions lor RUC md Eta forvcinu 
vs\ RUC analysis data for 40 analyst* cycles. 
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Figure 8: f\ valuation of Data Quality of the 48-ffour 
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Figure 9: Summary of tcnijic ratlin* differences for 
HOC and Eta forecasts to. RVC analysis (*F). 
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Figure 10: Summary of viod strength difttrt’nm 
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We observe that, based on the St D, one could state that 
RUC forecasts are marginally better than the 
comparable Eta forecasts for the same duration into the 
future up to 12 hours, at which point they are the same. 
But one can also see that 6 hour Eta forecasts are as 
good or better than 12hour RUC forecasts. The Eta 
forecasts maintain almost the same relative quality up 
to the maximum duration of 48hours tested (and by 
extension to the 60-hour duration available from Eta). 
Observe that the StD of temperature differences stays 
below 2 degree Fahrenheit for forecasts up to 48 hours, 
and StD of wind strength differences stay below 1 lft/s 
(6 .5 knots) for even the maximum duration forecasts. 
The wind direction difference, although with mean 
about zero, has a StD of a steady 50 degrees. These and 
some more detailed observations form the basis for the 
proposed switching protocol when either RUC alone or 
both RUC and Eta data become unavailable 

CONCLUSIONS 

The purpose of this study was to verify the quality of 
the Eta forecasts and to propose a decision algorithm 
for switching from RUC to Eta in case of outages. First 
let us state that Eta forecast quality is not in doubt in 
view of the results presented in the previous paragraph. 
Next, the protocol for switching will be based on the 
presumed quality of the forecasts. For example. If 
access to both RUC and Eta Is denied simultaneously, 
one should continue running CTAS with RUC forecasts 
as long as possible (from 9 to 12 hours, depending 
when in the RUC cycle the outage occurred). If RUC is 
interrupted, but not Eta, then one should switch at the 
next Eta update that would be at most 6 hours after the 
interrupt There are some additional considerations: 
CTAS expects a new RUC weather file to be made 
available every hour. When interrupts occur it may be 
necessary to create hourly files by interpolating from 
two adjacent forecast files. This would be true for RUC 
after the third hour and is always true for Eta forecasts 
as is clear from Figure 2 showing the RUC and Eta 
output products. The complete algorithm depends on 
the exact time of start of the outage compared to the 
underlying RUC and Eta cycle time, but it is based on 
the presumed quality of RUC and Eta forecasts at any 
given time. 
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ABSTRACT 

This paper addresses the issue of developing weather- 
impacted routes for the Final Approach Spacing Tool 
(FAST). FAST relies on adaptation data that includes 
nominal terminal area routes and degrees of freedom to 
generate optimum landing sequences and runway 
assignments. However, during adverse weather some 
adapted routes may become unavailable due to the 
presence of hazardous weather. If FAST continues to 
generate trajectories using these routes, its schedule will 
not be accurate during the adverse weather. The 
objective of the study was to determine methods for 
incorporating severe weather products and weather- 
impacted route data into FAST. 

BSTRQPVaiQN 

Increasing air traffic demand in the face of limited 
runways and airspace has made improving the 
efficiency of the nation's air traffic system one of the 
Federal Aviation Administration's (FAA) top priorities 
[1]. New decision support tools are being developed to 
assist Air Traffic Controllers and Traffic Managers in 
achieving these efficiency gains while maintaining 
safety. 

One such tool is the Final Approach Spacing Tool 
(FAST), which is an element of the Center/TRACON 
Automation System (CTAS) being developed by NASA 
Ames Research Center. There are two versions of 
FAST: Passive FAST (pFAST) and Active FAST 
(aFAST). Passive FAST provides controllers with 
runway assignments for delay balancing and landing 
sequence numbers to optimize the landing order. 
Active FAST additionally provides controllers with 
heading, speed and altitude advisories to achieve these 
optimal sequences [2]. Passive FAST has been 
operationally tested at the Dallas/Ft Worth 


•This work was performed for the National Aeronautic* and Space 
Administration under Air Force Coo tract No. FI962S-CO-OOOQ2 
tCcpyright C 2001 by M.l.T. Published by the American Institute 
of Aeronautics and Astronautics, Inc., with penrmsioo. 


International Airport (DFW) and has demonstrated up 
to a 13% airport throughput increase [3]. 

Active FAST, which is currently under development, is 
expected to offer additional increases in airport 
capacity. For the remainder of this paper, "FAST" will 
be used gencrically to refer to both systems. 

FAST relies on adaptation data that includes nominal 
TRACON routes and degrees of freedom to generate 
optimum landing sequences and runway assignments. 
However, during adverse weather, some adapted routes 
may become unavailable due to the presence of 
hazardous weather. If FAST continues to generate 
trajectories using these routes, its schedule will not be 
accurate during the bad weather. The objective of the 
study was to determine methods for incorporating 
severe weather products and weather-impacted route 
data into FAST. 

Air Traffic Control (ATC) decision support tools have 
many sources of weather data. Table 1 outlines some 
of these data sources and the information supplied. 

Since FAST operates in the terminal area, the key 
source for weather information for this study was the 
integrated Terminal Weather System (ITWS). ITWS 
was developed under FAA support by MIT Lincoln 
Laboratory (MIT/LL) and is currently being 
implemented for deployment at 35 locations across the 
continental United States. ITWS integrates weather 
data from a variety of sources, including Terminal 
Doppler Weather Radar (TDWR), Next Generation 
Weather Radar (NEXRAD), Airport Surveillance Radar 
(ASR-9), Rapid Update Cycle (RUC) model and 
surface sensor information. It generates a variety of 
weather display products for Traffic Management 
Coordinators (TMCs), controllers, airlines and other 
users. These products have been operationally tested at 
several sites including Dallas/Ft. Worth, Atlanta, 
Denver, Los Angeles and Miami, and have 
demonstrated significant operational benefits by 
allowing FAA personnel to improve their ability to 
manage traffic during bad weather [4]. 


57 



Ai-Vt Guidance, Navigation <k CofUroJ Conference Mootrai, Quebec, August 6-9, 2001 


T h»5 study focuses on the integration uod use (*>t‘ various 
ITW b weather pioiiucis in FAST, it alto addresses the: 
dbpiny of an integrated system to the end-user. ijj &hu: 
ense. tdcittiffed use?;; of tfrc system would «>e the Traffic 
Manager* m either tlvc TftACON 1 ot foe Air Route 
Traffic Control Center (A.RTOO. The italicized fTW'S 
product.? from Tahl-t I vrf*l b« cmpjv*?i2jc:i m tins proper. 

mBQAQl 

A study by Rhode; and Povvlak f S] revealed pilot 
febhavior m choosing to ^rtetrate or devote around 
to&ixioas wssther depends on several factors. 
Conventional wisdom is that piiotx tteviute around 
Video iotegrasor atKl Processor (Vi?) level 3 or high ex 
wetter. Howe vet, ixt evaiarnmg 1 ,052 encounters with 
v^afocr a- ;>FW, ft was found tl'at a significant 
sroponfoo of the encounter*; (>Ui%3 reunited in 
pemstration of VIP level ]% 4 or 5 weather, It was also 
formcl that the closer &» aircraft was to the airport, she 
mors likeiy it would penetrate VIP lave) 3 or higher 
weather A neural net classifier was developed to 
create Probability of Deviation {RODIr.V) maps giving 
die probability of a typical aircraft penetrauny. a given 
ofgioti of airspace. 

vVuh the development of the Tex rfo reft Convncftvu 
Weather forecast ftTCWPh Woifson, so at. }‘6‘] showed 
a significant increase in sk?ll could toe achieved m 
tweea-iting the location of V|p level 3 and higher 
weather for fine storms up to one hour in the future. 


TC.WF has been rumiing txperimeateHy at DFW since 
the summer of !9')V. This product is generated in 
conjunctior* with the ITW8 tcUbed operated by 
Ml VLL at DFW. The TO WE ha>. been pixipiwed as a 
Prs-ftbnned Prcxtect Improvement (V31) to 17 WS 

By applying the TCWF elgoritirtn to the forecasted 
weather products irt ITWS, tip to one hour of predicted 
VIP level 3 or higher wemhtr probability man* at can 
be generated at IsVmmnte ftitetvab. The neural net 
classifier used to create the i ; OD£V maps can fltem be 
applied to the TCWF forecasts to generate forecasted 
PODPV (FPODEV) maps. The PODEV end PRODEV 
maps arc ovortayed w;th the PAST adapted aircraft 
routes; to generate the probability of weather impact for 
each rnu-c segment from die present t : > one hour in the 
future in 1 0-minute inlet vah. 

An experimental weather s mi traffic Graphical User 
Interface (GUI), referred so as the Offline Traffic and 
Weather Display iOTWD), was developed to facilitate 
use of the weather* inspected route information. The 
aircraft, weather products and FAST aikptation touting 
information are displayed on OTWD. The roue 
segments -uc oolof-coded to reiteM the probability of 
whither impact at a pari ten tar time. A time slider was 
created to allow- the display of the w&lhur forecast jl \ id 
show the weather impact on aircraft route 
segments. Tile met ono also select any route segment to 
determine the onset and deration of hazardous- weather 
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An experimental version of the Aiiport Configuration 
Manager (ACM) was developed to facilitate TMC use 
of the weather-impacted routes information. ACM 
shows the impact upon the schedule of weather- 
impacted routes by predicting future airport 
configuration changes for each adapted aiiport This 
capability alerts the TMC to weather changes in the 
terminal area that could affect the operation of the 
airport up to one hour in advance. It also suggests an 
appropriate future airport configuration based on 
predicted weather data. 

Figure 1 demonstrates possible uses of ITWS weather 
products in FAST, from acquisition and parsing of the 
data to predicting impact on aircraft and terminal area 
operations. ITWS created precipitation maps can be 
converted into PODEV maps, which then can be used 
to identify routes and aircraft that are impacted. Delays 
caused by the hazardous weather can then be 
considered for ATC decisions such as airspace closures 
and re-routes. In addition, severe weather forecasts can 
be used to predict impacted aircraft and routes up to an 
hour in the future, as well as help determine appropriate 
future airport configurations. Schedule changes based 
on the weather forecasts can be utilized by TMCs to 
facilitate further ATC decisions. 


RESULTS 

OFF-LINE TRAFFIC AND WEATHER DISPLAY 

The identification and display of weather-impacted 
routes has been developed and demonstrated both on 
OTWD and on the FAST Plan view Display GUI 
(PGU1), sec Figure 2 [3]. Initial evaluation was 
completed using feedback from MIT/LL in-house 
controllers resulting in die current GUI design. The 
development of OTWD promotes the operational 
concept evaluation by allowing for the rapid 
prototyping of display modifications. However, the 
current functionality has only been implemented as a 
playback system. Additional work is needed to 
implement these functions into a real-time tool and to 
evaluate the initial operational concept 

OTWD displays aircraft and FAST terminal area 
routing, as well as the various weather products. The 
tool reads FAST generated playback files and 
synchronizes its internal clock to die data file time. As 
the weather products are read into the system at the 
appropriate time, OTWD analyzes the routing structure 
and weather data to determine the time, duration and 
severity of impaction on each route segment These 
route segments make up the nominal paths aircraft use 
to fly from the meter fix to each runway threshold. 


Probability of Current 

Deviation Impacted 

(PODEV) Map Routes 



Predicted Predicted 

Routes Schedules 


4 H M H ► 

iTK'SP’J CTAS WxFGUI FAST ctas future 


Figu re I. Weather Impacted Route* Proceislng 
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Figure 2 shows an usampte oi a PODEV nop for ibc 
DFW TRACX3N at 1 5.5'i GMT. The contour* range in 
intensity levol from 20*& • 100% shown by mortaring 
intensity of shaded area;* Notioc die selected touts 
Matting at die fk "BAMBB”. 1i bus a 
predicted duration of impact from "Current" to I6;0). 
with a Probability <;t Deviation of 80%. This means 
ihut fron 15-53 GMT imlii 16:0* GMT. a typical flight 
could be expected to deviate from this route segment 
550 % of the time. The information displayed for the. 
route segment beginning at “OO VIE** ahow$ tluit 
suiting ut 1 6:0 1 GMT tinough 16:51. GMT it w?U have 
a PODLV of 80%, though this route segment is not 
currently impacted by severe weather. 

1 he future weather forecast used to predict tfte airspace 
PODEYs can also be displayed on OTWD using the: 
slider bar, as seen in Figure 3. This figure shows the 
PPODEV at 16:14 GMT with tlx traffic data thorn 
15:54 GMT alkiwitig the user to visualize bow airspace 


problems might develop over ;.i me Though not 

currently pari of OTWD. tlx display oi the predicted 
aircraft position*; ai die "iooh ahead 0 lime has been 
denximtmtcd and could be beneficial 

Since PODEV maps are derived from 6 VIP level a nd 
TCWP data, OTWD b also capable <>t displaying these 
data. Figure 4 shows the TCW P tor the DFW 
TRACON The shaded areas correspond !<* predicted 
areas of VIP level i or higher weather activity twenty 
nntwtc* hon- tlx current time 

One area of research is «:• extend TCWF 10 other areas 
of tlx country. TCWP has Ixcn tected imd work:* well 
with Convective stumm. However. some area*? of the 
country tend k> have more localized thunderstorm cells. 
Forecasting these tyjxs of weather cell* is still being 
researched. Hence, die PODEV and terminal art*i 
forecasts for some areas of the country may run be 
mature cnouvh to predict impacted routes accurately. 
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AIRPORT CO WICKI^TIQti MMiACSB (A£ M i 

Using the identified impacted routes, ar. algorithm to 
help identify the most appropriate airport configuration 
for the turrom and predicted weatlier scenarios wai; 
developed, FAST ix&cik the concept of mr\x>n 
couftgataocm not only i<> identify ilie actual tun way » 
being nacd for landing* but a to to define the different 
procedure* used to control airers ft through the 
TRACON. Taking die current anti foixarasted weather 
mfortnatioi) into account, ACM ik designed to identify 
appropriate tiituic configuration changes urnl determine 
their time and duration. 

TIkt mos i important aspect fur dewrouning the 
appropriate airport configuration is the condition at the 
runway threshold. Thu includes the strength and 
direction of the wind. However, the H'WS weather 
data does not currently include terminal area wind 
prBiliotionv Therefore, Ihc search tint future auport 
configuration) is limited to the cimtieumiksn* adapted 
for FAST with runways landing in Ihe wane direct von as. 


the current configuration, The algorithm t«cessati«y 
assume* that the wind direction will remain relatively 
constant. 

Pint, the current airport configuration is identified. The 
adapted ixuifiguraiiottii aic seaiciied for others that 
contain any of tire runways u-sed in the current 
configuration. Thid forms our basic list of choice) tor 
future configuration changes. Ait runways used by any 
configuration in this list comprise our selection *>f 
possible future rurvwuys 

Next, the time when the current runways in use will 
change needs to be determihod. To do so, the adapted 
final approach aegmenc for each runway from our 
possible future runways iUt is examine*! For each 10- 
minute forecast interval over the next hour, it is 
determined whether each runway in urq>acted ot no*. 
Once the list of forecasted ik^ impacted runways, 
differ* from the runway* utilized m the current airport 
configuration, a future cuofiguraiiori change « 
identified 
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In practice the changing of airport coofign«*oon* cat) be 
difficult to perform. Ax\ arrival aircraft \u assigned to a 
runway upon catering the TRACON. An aircraft can 
take up to 20 minuter to tty through ti>c TRACON 
under normal condition*, so changing the an port 
configuration loo otkn would not be practical. 
Therefore future configuration changes are not 
automuted n> the ;>yMem. Instead a possible change is 
indicated, and the Traffic Managers can use the 
information to determine the bent course of actum. 


ITWS wind data to augment the RUC force**! data. The 
inclusion of additional weather products would be 
siraightfotwaxd to implement. RGDKV maps arc 
cun cully not included in ilk* operational ITWS system 
arul either another source of these map* would need to 
be created, or the coftwaie to create the K>DEV maps 
would need to be integrated into FAST. Once the 
impacted route software is m place, the value of AMC 
could be evaluated using operational TMCa at the f )FW 
TRACON 


NEXT STEPS 


Immediate next steps include implementing u reuUhne 
version of the OTWi) ut DFW. One Component of the 
implementation would be to create the infrastructure 
and interface to the i cal- time ITWS data at DRW 
Integration work has already begun for the ingestion of 


*1 lie use of ITWS wind forecasts to help manage the 
airport configuration con id also be rmtartxg*ied As 
stated above* ACM use* die tunways available for 
uircrat* landing in tire curreni direction to determine 
future runway uwage. However, it would be more 
uaeftil to include prediction;) of changes m the wind 
dirceti'W (hat fxuld impact th<? direction aircraft can 
l*nd. 
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As an extension to ACM, the identification and 
notification of impacted meter fixes to FAST could be 
added. This would allow the Traffic Managers to close 
meter fixes up to an hour in advance, giving aircraft in 
the En Route airepace the opportunity to re-route. 

OTWD could also be expanded from the terminal area 
to the en route area. With the incorporation of other 
weather data sources such as turbulence maps and cn 
route weather forecasts, the same tools could be applied 
to help predict airspace closure and aircraft intent 
outside the terminal airspace. This could improve 
predictions necessary for conflict detection and airspace 
overloading. 

Once the prediction of weather-impacted airspace and 
aircraft routes is complete, a logical next step would be 
to identify specific aircraft with impacted trajectories. 
This activity would tie together the predicted 
configurations and weather-impacted routes 
development as it pertains to the runway assignment of 
aircraft in the terminal area. Aircraft could be color 
coded by weather-impact to alert the controllers to the 
need to utilize alternate routes. 

An initial implementation could allow TMCs to 
graphically modify weather-impacted routes to create 
alternate nominal route segments for FAST to use for 
determining the path of flight in the terminal area. This 
would side step the procedural complexity of managing 
arrival and departure airspace. 

Based on route-building algorithms and heuristics 
investigated by Krozel et. al.[7], the indirect effects of 
severe weather on the closure of airspace due to 
thunderstorm downstream turbulence and severe 
weather in missed approach areas could be researched 
[8]. From these analyses, FAST could ultimately 
develop the alternate routes automatically. 

CONCLISIONS 

A method for connecting current and forecasted 
hazardous weather with FAST has been developed. 
The method utilizes previous work for creating 
probabilistic maps of pilot behavtv, * i penetrating or 
deviating around hazardous weal* This work has 
been extended by applying a new method for 
forecasting up to one hour in advance the probability of 
significant weather in a given region. These two 
methods are combined to produce forecast maps of 
probability of deviation for up to one hour in advance. 
These probability maps are applied to FAST adaptation 
data to determine current and forecasted weather- 
impacted routes. These weather-impacted routes are 
used in an intrractive display that could be used to 


provide guidance to Traffic Managers in choosing 

airport configurations in the presence of hazardous 

weather. 
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USING SURF AC E SURVEILLANCE TO HELP REDUCE TAXI DELAYS* 


Jorry i> Welch and Steven R. DussoUtri 
KUT Lincoln Laboratory 
244 Wood Strmtl 
Uxmgtork MA 02420-9100 

ABSTRACT 

Taxi delay is the largest of a!) aviation movement 
delays. However, taxi-out delay* tiave not received 
attention equal to that focused on airborne ddaya 
because teXMMit delays often result from downstream 
problems. Also, until recently, ilicrc was no practical 
mean* of tracking surface movements. New surface 
surveillance technology will revolutionize surface 
management by providing data, for planning, liming, 
buchJ roomtoriag surface operations. This (raper |>r»[>ose:j 
a simple aid to help managd departure taxi queues ami 
help exploit existing departure capacity, while avoiding 
the delays that resuit from saturated queues and 
unbalanced runways, The proposed decision aide will 
use archived surveillance data to quantify queuing 
behavior and model departure capacity, and it will use 
retd 'time surveillance to truck capacity changes and 
monitor d*e state of the taxi queues. 

Ta> i delay result* in the hugest direct operating cost to 
US uk carrier# of all delays. Fig. I shows that the 
average taxi out delay in minutes per flight is 
approximately twice the airborne delay/ Although 
aircraft bum fuel roughly .5 times tksier when abbot ne, 
crew and equipment costs make the upend rate for 
taxiing aircraft about 2/3 vital for airborne aircraft. 
Consequently, the coat of taxi out delay exceeds that 
for airborne delay by about 1 / 3 . totaling more than one 
billion dollar* amnudty. 

Oti avrra gt\ tax; -out delay is 3 nines large r then taxi- in 
delay. This ailuatlon suggests that surface aidu vrilt 
likely focus first on depan ultra. Mon! taxi oul delays are 
associated with surface queuing ppx.e\*es that arc 
visible lo the curiae? surveillance system, in contra?.*, 
matiagefoetd of the taxi-in process would require j*cccs* 
to airborne surveillance data. Therefore sotuc* 
planning ai d* for arrival * would need to be integrated 
that focus •>« with ARTS terminal surveillance data or 
the Center TRACflN Automatic) System (OTAS;. n 

Mhu % r<y\ >'*'1 pcjfcr^K-i for the ?<**fcn«} AcrooaalKA :w*i .V*»* 
Xdtramjunf W tauter Aif T-V**.* Coci'.*k! Ftr>s2MX)*0 •*#*»? 

• fVsp; rvhr 0 2fi0j by M l T. Pubhihca by the. An\ers*n ht>»iu*!r* ••{ 
.V*r>#;vi*tv . A*tr;f*v;*:ir., fisc . riffc 


Stepitert C. Atkin* 

VA£4 Ames Research Center 
M alt Stop 210* 

Mqffeit FkLi> C d 9403**1060 

Surface aids departures could rely ooteiy or* surface 
surveillance *n<f flight plan data, and could operate 
independently of ARTS mod CTA5. 



TAXI our TAX? *M 


f igtwt? L Tbe cost of delay - U«£, airline fleet cost 
estimates for 1999 from { A VA> 2000)/ 

Efforts to xeduce taxvout delay have not been as 
vigorous as those f<>cuiaxt on airborne delay . One likely 
reason is iha* taxi out delay often stems from 
downstream problem? such as capacity limitations in 
terminal and en route airspace as well nt> at Jbc 
declination airport. Delay stylistic* show a 10% 
correlation between (as i out delay for depatuites trom 
DRV and airborne delay for flights to DFW/ This 
mean* that, on day a when arrivals incur longer delays, 
departures at the came airport nsualiy take more time to 
get off the ground. Consequently, there has been a 
sense dim I idle can be done <m the surface to reduce 
mx> out cost. 

The airport surface has been ihc only dome&ic aviation 
domain without an automatic means of digitally 
tracking and identify ing aircraft. There ha-; tieou no 
practical mean? of tracking surface movements to 
uiid<:rs:at:d or contra; the taxi process, oiJUcr than luang 
reported push tare* and wheel* off tutaess. 

Hic : h petftxmanct cooperative surface surveillance with 
mrcrafl identification has now been demonstrated/** 
This new technology tvill provide an opportunity to 
revolutionize air traffic manage men! on the airport 
surface. There, has been considerable research in the 
application of surface surveillance tor enhancing 
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surface safety.* 19 There has been some research in 
analyzing and understanding taxi delays. 11 An attempt 
was made to mitigate taxi delays by predicting queue 
lengths from schedule data. 15 However, no work has 
been reported on the use of surface surveillance to help 
enhance surface efficiency and reduce taxi delay. 
Surface surveillance will help improve the 
understanding of taxi delay mechanisms and will 
provide data for characterizing operational constraints, 
planning taxi operations, timing taxi clearances, and 
monitoring aircraft compliance. 

The airport surface has characteristics that make it 
attractive for implementing automatic surveillance- 
based decision support tools. The surface is the only 
aviation domain in which most aircraft follow a 
relatively small number of rigidly constrained paths. 
This may make it feasible to use surveillance 
information alone to evolve a database that fully 
characterizes the geometry of the taxi paths. Generic 
adaptation algorithms could be designed to use this taxi 
path data to automatically adapt surface decision 
support tools to each new' airport site. Automatic 
adaptation would significantly reduce the cost of 
implementing surface decision support tools relative to 
terminal and en route tools. 

It appears feasible to develop operational software that 
uses surface surveillance information alone to track 
runway configuration changes in real time. By adding 
flight plan data to the surveillance data, automatic 
algorithms should be able to determine the current 
surface queuing delay status for each runway and 
aircraft. Other databases could relate taxi delay to 
demand and relate demand to the day of the week and 
the time of day. Operational algorithms could use these 
delay and demand databases in conjunction with 
measurements of departure queuing status to 
automatically predict the near-term departure 
throughput for each runway. With the addition of 
surveillance data for arriving aircraft, operational 
algorithms could also predict near-term arrival 
throughput 

A well-designed surface surveillance system can 
provide complete coverage of the important delay- 
controlling queues on the airport surface. Surveillance 
coverage of airport surface queues can provide the 
visibility into surface traffic flow that is essential to 
close the loop on suggested control actions. Surface 
surveillance provides a means of determining departure 
throughput performance, h allows unambiguous 
determination of the departure runway of each aircraft 
By contrast, the ARTS surveillance system cannot 
reliably associate departures with runways because of 
its low-elevation coverage cut-off. This is particularly 


true when the departures turn immediately after 
takeoff. 5 

A simple taxi-out aid could use this surveillance 
information to predict the queue length for the next few 
minutes and advise optimum near-term target pushback 
rates. Such advisories would have the advantage that 
they would not require any manually generated 
information from controllers or aircraft operators. 
There would be no need for controllers to provide the 
currant runway configuration and no need for pushback 
predictions from aircraft operators. 

A SIMPLE TAXI-OUT AID 

The most elementary taxi aid would display surface 
traffic with aircraft identities to all parties interested in 
surface traffic management. Distributing surveillance 
information with aircraft tags would likely benefit 
aircraft operators and controllers and would require 
little research other than finding means to manage the 
tags to avoid display clutter. The functionality 
envisioned in this paper goes further and provides 
decision support aids based on surface surveillance that 
would attempt to directly help controllers and air 
carriers work together to improve surface movement 
efficiency. This is done by providing a simple 
visualization of the queuing situation along with 
pushback advisories to optimize the taxi-out process 
and balance runways. 

At airports with multiple runways, surveillance data 
alone presents a good tactical view of the airport, but 
sometimes paints a confusing picture of the strategic 
situation. Early in the taxi-out process, it is not always 
clear which runway each aircraft is heading for. 
Although the total number of taxiing aircraft may be 
apparent, it is difficult to see the queue lengths for 
individual runways. The surveillance display also does 
not provide information on predicted taxi times. 

Experience with the Surface Movement Advisor (SMA) 
program at Atlanta’s airport conclusions from current 
NASA-sponsored research on causes of departure delay 
and inputs from aircraft operators who routinely 
experience departure delays all suggest that an 
automation aid to help visualize and control taxi queue 
lengths for departing aircraft would benefit both 
controllers and aircraft operators. 3 ‘°* ,6,lx,4 * ,A5 * 9,a 

This paper proposes such a simple taxi-out aid to help 
with queue management. Figure 2 shows a notional 
operator interface for the aid when used at an airport 
with dual departure runways. The interface is presented 
merely as a concrete illustration of the proposed 
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udVmmtion water*. it has not beta prototyped or 
trswd fo r operator acceptability. 

ft display* the slat us at the laxiway system u.k! depicts 
the queuing situation at multi-runway airports by 
associating each aircraft with a runway and ifcpicting 
the airport \n a simplified map format similar to those 
*.»cd fo** mapping subway line*. It *how* the loading of 
the taxi paths as welt as the runway queues, indicating 
the progress at all aircraft taxiing towards tuck 
departure runway. 

Ihe runway maps are divided ramtinaliy into Pushback, 
Taxi, Runway Crowing Queue, and Runway Queue 
notions. The taxi and Queue sections contain bias 
consisting of columns trf boxes indicating numbers of 
aircraft Bach box corresponds to an dreraii. Clicking 
<m a box displays a dau* nig with the aircraft ID. In t he 
figure, tun way I8C has two planes waiting to take off 
In its runway queue. Where a taxi path crosses an 
active runway, the graph ia split by a crc*Bittg queue sw 
shown ixrr Runway l#R. 

The Taxi section la shown with M> him;. Typically, 
each taxi inn contains rso more than one aircraft. Large 
airports use more taxi bins to keep die individual taxi 
bia occupancy to a single aireratr. 


A Crossing Queue section has a single cotonro showing 
die number of aircraft waiting to cross the active 
runway. When an aiutaft first tmets a runway crossing 
queue or x runway queue, its box changes color to 
Jndicak ?» it is now under local corned Finally. the 
graph for each rtmtv*y er«h in a Runway Queue section 
containing a single column showing the numbet of 
aircraft waiting at the runway entrance for takeoff 
clearances. 

Bach runway also has a pushback rate advisor. When 
the departure demand is low, there is no need to limit 
the pushback rate. Thud, although the display operates 
continually, its principal benefit occurs ix« departure 
push conditions. T he height of the b «r recommends a 
pushback rate for that runway. A positive bar indicates 
that the pushback rate can be increased, and a negative 
i»r advises a reduced mte. hi a departure push the gc&l 
is to fxro the bars for both runways. Zeroing the bars 
optimize* the pushback rates and equalises the tax! 
times for ihe last aircraft in eac h of ibe runway queue-;. 
(Note that -he las- aircraft in a runway queue may not 
be the most recent aircraft to push back for that 
runway.) The display also indicates the estimated taxi 
time for tiw last aircraft in each runway queue. 


PusiUw bar advises tac.re*is«rf maftbork <«u 



Sbn w* numb* or Dianes at am way entrance i 


18 R 


raKl quuuo may be ipfit hy CrosDng Queue?*) 








Taxi 


26 . 


IT push beck tuo tost, b« nous neyatfye, 

Og* t* to mo the twir. 


Figure 2. Nominal display for m» airport with dual departure runways 


67 


p 

A1AA Guidance, Navigatioc & Control Conference, Montreal, Quebec, August 6-9, 2001 


Equalizing the taxi times for the last aircraft in the 
queues balances the runways despite differences in their 
overall taxi path lengths or differences in their 
departure capacities. Runway capacity differences can 
be caused by departure airspace constraints or by 
differences in mode of operation. For example, a 
runway dedicated to departures has greater departure 
capacity than one shared by arrivals. Balancing die taxi 
times assures that the flow to one runway does not dry 
up earlier than the flow to the other at the end of a 
departure rush. The simple taxi-out aid does not 
attempt to balance taxi delays. Balancing taxi delays 
(where delay for each runway is reckoned relative to its 
minimum taxi-out time) could assign too many aircraft 
to the runway with the longer taxi path. 

USING THE TAXI AID 

Large hubbing airports with multiple departure runways 
often allow airline gate/ramp control personnel to 
control the pushback and ramp taxi process for their 
own aircraft They interact with the FAA Gearance 
Control. Ground Control, and Local Control positions. 

After FAA Gearance Control has issued a flight plan 
clearance and handed off the flight strip to FAA 
Ground Control, the next step in the departure process 
occurs when the pilot notifies the carrier's Gate/Ramp 
Control that he is ready to push back. Gate/Ramp 
Control issues the pushback clearance and clears the 
aircraft to the desired taxiway entrance spot Upon 
reaching the spot, the pilot monitors the Ground 
Control frequency in anticipation of a taxi clearance. 

Ground Control positively identifies each aircraft that is 
ready to enter the surface movement area and, when 
ready, clears it to enter the taxiway. Ground Control 
handles the aircraft until it first approaches an active 
runway. Transfer to FAA Local Control may take place 
before crossing an active runway or before entering the 
departure runway, whichever occurs first Ground 
Control hands off the flight strip to Local Control and 
authorizes the pilot to switch to the Local Control 
frequency. Local Control then issues all remaining 
clearances including the takeoff clearance. 

If two or more carriers simultaneously push back 
aircraft for departure, FAA Ground Control establishes 
the overall taxi-out order to sequence the aircraft on the 
runways in a fair and equitable manner subject to 
departure flow control and wake vortex spacing 
restrictions. Ground Control also exercises short-term 
control over the runway balance and the queue lengths 
for the individual runways by managing taxi out from 
the ramp area. 


Surface surveillance displays will be used in the tower, 
the Traffic Management Units in the TRACON, the en 
route Center (where it could help predict near-term 
departure demand), and the airline operational centers 
and gate and ramp control positions. The simple 
taxi-out aid will be used at all the same locations. 

In departure pushes FAA tower controllers estimate the 
size of departure queues by monitoring the occupancy 
of the flight strip bays and by visually observing the 
aircraft on the runways and taxi ways. The taxi-out aid 
provides Ground Control with a direct estimate of the 
current runway balance and suggests the optimum 
apportionment of aircraft to multiple departure 
runways. By comparing the queue size to the visual 
scene or to the tower flight strip bays its users can 
assure that it accounts for all the aircraft, assigns each 
aircraft to its proper departure runway, and accurately 
shows the distribution of the aircraft along the paths to 
the runways. Users can dick on aircraft boxes to 
display data tags confirming aircraft identity and status. 

Gate/Ramp controllers often have difficulty directly 
observing the details of the other carriers' departure 
operations or the overall departure demand and 
capacity. The taxi-out aid provides the needed 
visibility and suggests near-term limits to the pushback 
rate for all of the carriers. Individual carriers can then 
infer their own pushback rate limits by comparing the 
pushback advisories with their recent departure flow 
performance. 

The simple taxi-out aid facilitates collaborative decision 
making among the carriers and ATC by providing all 
participants with common situational awareness. They 
collaboratively determine the throughput performance 
objectives for the airport and enter the desired departure 
rate into the simple taxi-out aid. They can change the 
control performance goals at any time, resulting in 
immediate changes to the displayed pushback rate 
advisories. When a single carrier dominates the 
departure push the pushback advisories more directly 
apply to the dominant carrier, simplifying the 
cooperative management of gate or ramp holds between 
the hubbing carrier and FAA Ground Control. 

TOE PUSHBACK ADVISOR Y ALGORITHM 

The taxi-out aid automatically determines the 
relationship between departure throughput, taxi-out 
time, and effective queue length for each 
taxiway/runway path. The effective queue for a runway 
includes all aircraft on their way to the runway as well 
as the aircraft in the physical queue that the runway 
entrance. The relationship between these performance 
measures is obtained by continually measuring. 
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tracking, and archiving each of these quantities. The 
taxi aid analyzes the archived data under all operating 
conditions to determine and model the needed 
relationships. 

Figure 3 illustrates the relationships used to develop the 
departure performance database. It relates the overall 
mean departure throughput to the effective queue length 
for BOS in August 1991, which includes all aircraft that 
have pushed back but not yet taken off. 



0 5 10 15 20 25 30 

Estimated taxi queue length at start of period 


Figure 3. Mean departure rate vs. taxi occupancy 
for BOS, Ang 1991 - from (Shumsky 1997). 1 * 

Taxi-out delays exhibit classical queuing behavior. In 
departure rushes, aircraft fill the taxiways, queues grow 
at runways, the ambient runway/airspace departure 
capacity limits the takeoff rate, and delay increases 
faster than the queue length. The service rate of a 
queue is always the minimum of the demand and the 
capacity. Figure 3 shows that when there art less than 
15 departure aircraft on the taxiways, there are gaps in 
the flow to the runway, and the takeoff rate is 
determined mainly by the departure demand. When 
there are many more aircraft on the taxiways than 15, 
there is usually a queue of aircraft at the runway 
entrances. The takeoff rate is then determined by the 
departure capacity of the airport The variability in 
mean takeoff rate for large queue lengths reflects the 
variation in airport capacity that occurred during the 1- 
roonth data-gathering period. If the operator desires a 
total throughput of seven aircraft in 10 minutes, this 
curve tells him that the queue length should be about 
fifteen. Increasing the queue length above 15 will have 
little operational effect other than to increase taxi-out 
delays. 

The taxi aid bases the pushback limit calculation on the 
principle that the rate of growth of the departure queue 
is determined by the difference between the pushback 
rate and the runway take-off rate. The taxi aid 
determines the pushback limit from the queuing 
conservation relationship: 


P-Nt-No + D, 

where P is the number of pushbacks during the next T 
minutes; N T Is the desired queue length T minutes in the 
future; N 0 is the present queue length, and D is the 
number of takeoffs expected during the next T minutes. 

The desired future queue length, N T is obtained from 
the performance database after a throughput goal has 
been established by the users. The present queue 
length, N 0 is obtained from surface surveillance. The 
number of takeoffs during the next T minutes is 
obtained by predicting the departure rate over the next 
T minutes. 

To predict the departure rate the simple departure aid 
uses historical data relating queue length to departure 
rate. The departure rate is determined by the departure 
capacity of the airport when the departure queue is 
large. In a departure push the departure capacity varies 
slowly and can be reliably tracked and predicted in real 
time with no knowledge of the weather, the runway 
configuration, or the arrival/departure mix. 19 Shumsky 
was able to forecast the takeoff rate 30 minutes into the 
future with a RMS error of about two departures per 
hour. He achieved his best capacity estimates by fitting 
an analytical approximation to data of the type shown 
In Fig. 3. He used an exponential curve fit to provide a 
static model of the relationship between queue length 
and departure rate. He then added a dynamic term to 
adjust the curve up or down at 10-minute intervals. The 
dynamic adjustment linearly tracked the observed 
takeoff rate by smoothing over the residuals of the 
takeoff estimate at each update interv al. 

Shumsky did not have surface surveillance data to help 
refine his estimate. Knowing the runway configuration, 
the distribution of aircraft on the surface, and the 
balance between arrivals and departures for each 
runway makes it possible to improve the static 
estimates by modeling families of curves for different 
operating conditions. Tracking repeatable daily or 
hourly variations can further improve the estimation 
process. 

The taxi aid also estimates the total taxi-out time for the 
next pushback to help balance the queues. The 
relationship between effective queue length and taxi 
time is also obtained from the performance database. 
This is illustrated in Figure 4, which relates the mean 
taxi-out time to the queue length for the entire departure 
taxiway System at Boston (BOS) from January to 
March 1997. This curve indicates that a queue length 
of 15 aircraft produces an expected taxi-out time of 36 
minutes. The mean taxi-out delay relative to the 
unimpeded delay can also be obtained from this figure: 
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the taxi delay is approximately 36 - 12 = 24 minutes for 
a queue length of 15 aircraft. 




0 4 8 12 16 

Departure tawway occupancy (aircraft) 


Figure 4. Dependence of mean taxi out time on the 
aircraft co out is the BOS departure taxiway system 
at the time of pushback - from (Idris et al, 1999). 12 

To estimate the departure queue length in real time for a 
multi-runway airport, the taxi aid associates each 
aircraft with a departure runway and automatically 
detects changes in departure runway assignment At 
pushback the taxi aid estimates the departure runv'ay 
for each aircraft from the filed flight plan for the 
aircraft and the current runway configuration. The taxi 
aid checks the position of the aircraft at each 
subsequent update against a site-adapted list of 
departure taxi routes for that configuration to confirm 
that the aircraft is still taxiing towards the same runway. 


Although the simple taxi-out aid needs only flight plan 
data and surface surveillance data to function, 
additional data types would be valuable for surface 
movement predictions. Pre-departure pushback status 
information and automatic pushback notices from all 
aircraft operators at the airport, final approach 
surveillance data from ARTS, CTAS data, and data on 
operational constraints in departure airspace would all 
improve the performance of taxi aids. 

These data sources would help extend departure 
demand predictions farther into the future to better 
support strategic planning. Pre-departure status 
information from aircraft operators would be very 
useful if it were reliably available for all departures. 
Taxi aids will require timely automatic pushback daft at 
an airport whose surface surveillance does not provide 
reliable coverage in the ramp area. It is essential to 
measuring the gate and ramp performance of the 
surface surveillance system as the first step in 
implementing a taxi aid at an airport. 


CTAS data includes arrival planning information as 
well as ARTS approach surveillance data. These would 
both be necessary for more advanced taxi aids that 
manage arrivals as well as departures. 

Other information would be useful in the future to 
refine the departure capacity estimates. Aircraft type 
information (which is included in the flight plan) would 
help determine wake vortex spacing requirements on 
take-off. 7 Although it is possible to automatically 
detect runway configuration changes in near real time 
from either ARTS surveillance data or surface 
surveillance data knowledge of planned configuration 
changes would obviously improve the predictability of 
departure capacity. 1 Finally, information on current and 
planned departure airspace restrictions would help 
improve both the capacity estimates and the archived 
data used for modeling constraints on departure 
sequencing and timing. 


The principal benefit expected from the simple taxi aid 
is a reduction In direct operating costs for aircraft 
operators from reduced taxi-out delay. The aid will 
also reduce environmental pollution and controller 
workload. It will achieve these benefits by helping 
controllers and aircraft operators determine when to use 
gate and ramp holds to avoid overloading departure 
queues at individual runways. At airports with two or 
more departure runways, it will help reduce departure 
delays by balancing the queues to prevent under- 
utilization of one runway while overloading another. 
Predicting such unbalances early will help reduce 
expensive delays for aircraft that are not yet committed 
to depart on a particular runway. 1 * However, when 
runways become unbalanced because of procedures that 
rigidly map departure fixes to runways, the solution 
will likely involve procedural changes to allow an 
aircraft to depart from a runway not normally 
associated with its planned departure fix. 

In addition to avoiding runway saturation and runway 
imbalances, it also helps ground controllers plan taxi 
clearances and assign departure runways in order to 
provide continuous streams of traffic for all departure 
runways. The magnitude of these benefits can be 
determined initially by analysis and then by using 
controUer-in-the-loop simulations and re-enacting 
demanding departure taxi scenarios. Baseline tests can 
be run on the scenarios to quantify the resulting 
departure delays with and without the help of pushback 
advisories. 

Although the initial functionality will be limited to 
pushback advisories and departure queue length 
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predictions, its use of surface surveillance and digital 
night plan data to predict surface movements will be a 
necessary first step in enabling more efficient 
management of surface movements and in supporting 
future functionality. Future efficiency -enhancing 
functionality based on the simple departure aid could: 

• help coordinate arrival and departure planning by 
accurately predicting departure times, 

• estimate taxi-in as weli as taxi-out delays for 
Individual aircraft 

• help provide runway assignment and sequencing 
advisories for arrivals as well as departures, 

• help support procedures that permit an aircraft to 
depart from a runway not normally associated with 
the departure fix designated in its flight plan, 

• help reduce runway crossing delay by reducing taxi 
time variances and excess buffers, 

• help improve reconfiguration planning, 

• help increase landing throughput by monitoring 
runway occupancy times, and 

• help reduce de-icing delays by monitoring de-icing 
queues and time violations. 

All of these efficiency enhancements will require the 
runway adaptation, surveillance and flight plan 
processing, surveillance data analysis, capacity 
prediction, and display generation capabilities that are 
required in a the simple taxi-out aid. These capabilities 
will also provide the basis for important surface safety 
functional ity. 

RISKS 

The fundamental risk inherent in any activity to develop 
decision support tools is the failure to identify the main 
operational problems and the consequential failure to 
deliver dollar, workload, or environmental benefits at 
the selected site. Data is needed to verify assumptions 
about delay and cost mechanisms at candidate sites. 
The surface surveillance system and FAA delay 
reporting systems must be used along with operational 
logs and consultation with controllers and aircraft 
operators to determine the principal causes of taxi-out 
delay. 

Failure to achieve user acceptance is a fundamental 
risk. The controller interface design and development 
process involves significant research risk. 
Familiarization, training, and display modifications will 
be required, and simulations will be needed with and 
without the pushback advisories and queue distribution 
graphics to determine if they are useful and reduce 
workload for the users. 


It will also be necessary to find space for a new display 
in the tower. Experiments with trial interface designs 
must begin as soon as possible in a tower environment 
The NASA tower cab simulator is ideal for this 
purpose. Adapting the NASA simulator to the chosen 
airport should begin as soon as the test site has been 
chosen. This adaptation is expensive, so there are 
budgetary risks involved in this site decision. However, 
evaluation of generic displays could begin immediately 
with a tower cab simulation of any airport 

Major research risks were noted above in discussing 
functionality and data interfaces. An important risk is 
the performance and accuracy of the key capacity 
tracking and prediction algorithms used to estimate the 
future departure queue length. These estimates will 
rely on processing and tracking routines that 
continuously monitor, record, and analyze operational 
surface surveillance at each airport Although data 
tracking is performed routinely in Air Traffic control 
surveillance systems, it is not commonly done to 
automatically obtain and update operational 
information needed for decision support automation . 3,10 
If it is not possible to predict departure capacity with 
sufficient accuracy solely from recent observations of 
departure queues and takeoff rates, it may be necessary 
to obtain additional sources of data. 

Important programmatic and technical risks involve 
access to, continued availability of, and coverage of the 
- surface surveillance data. Surveillance data is key to 
automatically determining the current queuing status for 
each departure runway and automatically predicting 
near-term departure demand. The FAA is developing a 
commercial surface multi lateratlon system intended to 
provide reliable surveillance and identification of all 
aircraft with operating transponders on airport 
movement areas. The data is fused with data from 
surface surveillance radars, with data from the ARTS 
terminal airborne surveillance system, and with flight 
plan data. 

The FAA has initiated development of an operationally 
deployable digital surface surveillance system known as 
ASDE-X. This program is intended to lead to an 
operational capability that includes multilateration in 
conjunction with a low-cost primary radar operating at 
X-band. 

A minor regulatory change is needed to obtain reliable 
surveillance and identification coverage for all aircraft 
in the ramp and movement a 'as. Transponder 
multilateration, which is the key to low-cost, clutter- 
free surveillance on the surface, depends on 
transmissions from aircraft transponders. Official FAA 
procedures must change at airports with multilateration 
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systems to mandate that aircraft on the airport ramp and 
movement areas leave their transponders operating at 
all times rather than turn them off as currently required. 
Because there appears to be no technical or interference 
problem from continuous operation of transponders on 
the surface of major airports there is no technical or 
operational impediment preventing this regulatory 
change. 30 

There is a small technical risk that surface surveillance 
may not provide reliable coverage in ramp and gate 
areas because of blockage of transponder transmissions 
by airport structures or by the airframes of other 
aircraft Recent tests at DFW addressed the coverage 
issue. 31 The gate coverage appeared good based on 
limited operation within the ramp and gate areas. 
Muldlateration coverage can always be improved by 
the addition of additional ground sensors. In addition, a 
few aircraft will not be visible because they have been 
intentionally wired to reduce controller workload by 
automatically switching off their transponders when 
they are on the ground. One of the early research 
activities in any program to develop surveillance-based 
taxi aids must be a complete characterization of the 
coverage issue. 
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LIST OF ACRONYMS 


ACM 

ARTCC 

ASR-9 

ATC 

Airpoit Configuration Manager 
Air Route Traffic Control Center 
Airport Surveillance Radar 
Air Traffic Control 

CONUS 

CTAS 

Continental United States 
Center/TRACON Automation System 

DFW 

Dallas/Ft. Worth International Airport 

ETA 

Estimated Time of Arrival 

FAA 

FAST 

FPODEV 

Federal Aviation Administration 
Final Approach Spacing Tool 
Forecasted Probability of Deviation 

GW 

Graphical User Interface 

rrws 

Integrated Terminal Weather System 

MTT/LL 

Massachusetts Institute of Technology Lincoln Laboratory 

NEXRAD 

Next Generation Weather Radar 

OTWD 

Offline Traffic and Weather Display 

P3I 

PGUI 

PODEV 

Pre-Planned Product Improvement 
Planview Display Graphical User Interface 
Probability of Deviation 

RUC 

Rapid Update Cycle 

TCWF 

TDWR 

TMC 

TRACON 

Terminal Convective Weather Forecast 
Terminal Doppler Weather Radar 
Traffic Management Coordinator 
Terminal Radar Approach Control 

VIP 

Video Integrator and Processor 
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